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Abstract: 

We introduce a duality between two-dimensional XY-spin models with symmetry-breaking pertur- 
bations and certain four-dimensional SU{2) and SU{2)/Z2 gauge theories, compactified on a small 
spatial circle M""^'^ x S"'^, and considered at temperatures near the deconfinement transition. In a Eu- 
clidean set up, the theory is defined on X T^. Similarly, thermal gauge theories of higher rank are 
dual to new families of "affine" XY-spin models with perturbations. For rank two, these are related to 
models used to describe the melting of a 2d crystal with a triangular lattice. The connection is made 
through a multi-component electric-magnetic Coulomb gas representation for both systems. Perturba- 
tions in the spin system map to topological defects in the gauge theory, such as monopolc-instantons 
or magnetic bions, and the vortices in the spin system map to the electrically charged W^-bosons in 
field theory (or vice versa, depending on the duality frame). The duality permits one to use the 
two-dimensional technology of spin systems to study the thermal deconfinement and discrete chiral 
transitions in four-dimensional SU{Nc) gauge theories with nf>l adjoint Weyl fermions. 
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1. Introduction, summary, and outline 

It is well-known since the late 70's that two-dimensional (2d) XY-spin models with appropri- 
ate symmetry-breaking perturbations map to a 2d Coulomb gas with electric and magnetic 
charges [1]. This beautiful duality permits analytic calculations of long distance correlation 
functions, phase diagrams, and critical indices in a large class of 2d statistical mechanics 
systems such as XY- and "clock"- (planar Potts) models (for reviews of methods and applica- 
tions see, for example, [2-4] and references therein). A vectorial version of the XY-model has 
been used in the study of melting of a 2d crystal with a triangular lattice [5]. The appropriate 
Coulomb gas is a system of electric (e) and magnetic (m) charges, and involves e-e and m-m 
Coulomb interactions, as well as e-m Aharonov-Bohm phase interactions [6]. 

In this paper, we introduce a long distance duality (equivalence) between 2d XY-spin 
models and certain four-dimensional (4d) gauge theories compactified on M2 X with 
boundary conditions as specified below. The connection is made through an electric-magnetic 
Coulomb gas representation for spin systems and gauge theories. As stated above, the con- 
nection of the spin systems with Coulomb gases is well understood. What is new is the 
realization that some four-dimensional gauge theories on M? X also admit an electric- 
magnetic Coulomb gas representation. This is the basis of the long-distance duality between 
XY-spin models and circle-compactified 4d gauge theories at finite temperature. 

Let us briefly summarize the main progress which makes the mapping of four-dimensional 
gauge theories to electric-magnetic Coulomb gases possible. The gauge theory we study is 
4d QCD with adjoint fermions, QCD(adj), formulated on X 

S^. Here, E>\ is a spatial 

(non-thermal) circle of circumference L and the fermions obey periodic boundary conditions. 
This theory does not undergo a center-symmetry changing phase transition as the radius of 
the is reduced [7] . This implies that the theory "abelianizes" and becomes weakly coupled 
and semi-classically calculable. Non-perturbative properties, which are difficult to study on 
M^'^, such as the generation of mass gap, confinement, and the realization of chiral symmetry 
can be studied analytically at small [8,9]. The non-perturbative long-distance dynam- 
ics of the theory is governed by magnetically charged topological molecules, the "magnetic 
bions" — molecular (or correlated) instanton events whose proliferation leads to mass gap and 
confinement. The vacuum of the theory is a dilute plasma of magnetic bions. The theory also 
possesses electrically charged particles, e.g., M^-bosons, which decouple from the long-distance 
dynamics at T = 0. The zero-temperature physics of QCD(adj) on X is reviewed in 
Section 2. 

We now turn to finite temperature, corresponding to compactifying the theory on ]R2 X 
S\ X S^. Electrically charged particles, e.g., VF-bosons, can now be excited according to 
their Boltzmann weight e~^"^^ . It turns out that their effect is non- negligible, as in studies 
of deconfinement in the 3d Polyakov model [10, 11] and in deformed Yang-Mills theory on 
M^'^ X S]^ [12]. Thus, near the deconfinement transition we must consider Coulomb gases of 
both electrically and magnetically charged excitations. This is the main rationale under the 
gauge theory /electric- magnetic Coulomb gas mapping on x S]^ x S^. This mapping reveals 
a rather rich structure that we have only begun to unravel. 

Below, in Sections 1.1, 1.2, and 1.3, we summarize the main results of this paper. 
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1.1 Deconfinement in SU(2) QCD(adj) 

We first review our results for an 5f7(2)/Z2 gauge theory with Uf massless adjoint Weyl 
fermions on X S^. We show that the physics near the deconfinement temperature is 
described by a classical 2d XY-spin model with a U{1) — t- Z4-breaking perturbation. This is 
the theory of angular "spin" variables 6x € (0, 27r], "living" on the sites x of a 2d lattice (of 
lattice spacing set to unity) with basis vectors fi, with nearest-neighbor spin-spin interactions. 
The partition function is Z = f T)9e~^^ ("/3ff" is used to denote the action, i.e. f3 below is 
not the inverse temperature), where: 

-PH= ^ |-cos(0,.+/i-0,) + ^ycos40, . (1.1) 

When matching to the gauge theory, the lattice spacing is of the order of the size L of the 
spatial circle. The spin-spin coupling n is normalized in a way useful for us later. 

We emphasize that the equivalence of (1.1) to the finite-temperature gauge theory is 
not simply an effective model for an order parameter based on Svetitsky-Yaffe universality 
[13]. Instead, the parameters of the lattice spin theory (1.1) can be precisely mapped to 
the microscopic parameters of the gauge theory, owing to the small-L calculability of the 
gauge dynamics. This map is worked out in Section 3, where the nature and role of the 
perturbative or non-perturbative objects driving the deconfinement transition is made quite 
explicit. We note that there have been earlier proposals and discussions of the role of various 
topological objects in the deconfinement transition, in the continuum and on the lattice, 
and that some bear resemblance to our discussion (see, for example, [14-21] and references 
therein). However, the study here stands out by being both analytic and under complete 
theoretical control. In addition, we study QCD(adj), a theory with massless fermions, and 
not pure Yang-Mills theory. 

In the lattice theory defined as in (1.1), the XY-model vortices map to electrically 
charged VF-bosons, while the Z4-preserving perturbation represents the magnetic bions of the 
SU{2)/'L2 QCD(adj) theory. The spin-spin coupling is expressed via the four-dimensional 
gauge coupling g4:{L), the size of the spatial circle L, and the temperature T, as: 

n = ^, (1.2) 

and determines the strength of the Coulomb interaction between the P^-bosons, while the 
(dual-) Coulomb interaction between magnetic bions is proportional to k^^. 

The global U{1) symmetry of the XY-model, Ox ^ Ox + c, is explicitly broken to Z4 
by the magnetic-bion induced potential term in (1.1). In terms of the symmetries of the 
microscopic SU{2)/Z2 gauge theory, the Z4 symmetry of the spin model (1.1) contains a 
discrete 7^2^ subgroup of the chiral symmetry of the gauge theory and the topological 
symmetry, which arises due to the nontrivial homotopy 7ri(S'f/(2)/Z2). The combination of 
these two symmetries accidentally enhances to give a I^J'^^ symmetry. 

There exists a lattice formulation different from (1.1) and appropriate to an SU{2) gauge 
theory (which has trivial vri, see the discussion in Section 3.3), where, instead of the topological 
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symmetry, one finds a center symmetry. In this case, the discrete symmetry group of 
the theory is x Tl^ ^ and does not enhance to a Z4. The reahzation of the symmetries and 
the behavior of the correlators of appropriate 't Hooft and Polyakov loops above and below 
the deconfinement transition are discussed in Section 3.3.3 and are summarized in eqns. (3.35) 
and (3.37). 

The lattice-spin model (1.1) is a member of a class of Zp-preserving models, defined as 
in (1.1), but with co's.ABx cospOx instead; these are sometimes also called the Zp "clock" 
models, because, in the limit of large y, the "spin" e*^ is forced to take one of p "clock" values. 
The Z4 model stands out in this class in that the critical renormalization group trajectory is 
under theoretical control, at small fugacities, along the entire renormalization group flow to 
the fixed point describing the theory at Tc- This makes it possible to obtain the analytic result 
for the divergence of the correlation length, given below in eqn. (1.3), see also Section 3.2 and 
Appendix B.2. The phase transition in the Z4 model occurs at k = 4, thus, from (1.2), the 
critical temperature is given by Tc ~ We note that K = 4 is also the critical coupling 

for the usual Berezinskii-Kosterlitz-Thouless (BKT) transition of the XY-model without the 
U{1) —7- Z4-breaking term. The transition in the Z4 model is also continuous, however, as 
opposed to BKT, it is of finite (but very large, at small fugacities) order. As we show in 
Section 3.2, the correlation length diverges as: 



as T —7- Tc from both sides of the transition. Here, yo and yo are exponentially small param- 
eters, essentially determined by the fugacities of the VF-bosons and magnetic bions at scales 
of order the lattice cutoff L~^. At small values of these fugacities, the critical correlators are 
essentially governed by a free field theory with BKT (k = 4) exponents. The phase transition 
in the Z4 "clock" model at k = 4 corresponds to the confinement/deconfinement transition 
in the SU {2){adj) theory. Both the string tension a and the dual string tension a vanish as 
the inverse correlation length (""^ as T —)• Tc from below (a) and above (a). 

An important property of the rank-one case is that the electric-magnetic Coulomb gas 
dual to the spin model (1.1) exhibits electric-magnetic (e-m) duality. This duality is not 
manifest in eqn. (1.1), but is evident from the Coulomb gas representation, see Section 3.3. 
It involves exchanging the fugacities of bions and I^-bosons and an inversion of the coupling 



Thus, the critical temperature Tc is precisely determined by the strength of the interaction 
K, at the point where the Coulomb gas is self-dual. This e-m duality property is shared by 
all Zp "clock" models. As usual with Kramers- Wannier-type dualities, it helps establish a 
candidate critical temperature. Using 2d CFT techniques, it has been shown [4] that, indeed, 
at the self-dual point the Zp=2,3,4 models map to known conformal field theories. For p = A, 
this is a free massless scalar field, even at large fugacities.^ 

^For completeness, we note that the p > 4 models have an intermediate massless phase, see, e.g., [22]. 




(1.3) 



(1.2): 



16 



(1.4) 



K 
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Phases of SU(2) and SU(2)/Z2 QCD(adj) 




L 



Figure 1: Symmetry realizations for SU{2) theory on X §i X S^: The SU{2) theory 
possesses x Z"^^ x SUirif) center, discrete, and continuous chiral symmetries (xS). For brevity, in 
each phase we only list the broken symmetries. A: Discrete and continuous broken. B: Discrete 
xS broken. C: Center symmetry broken. D: All symmetries broken. Symmetry realizations 
for SU{2)/'L2. theory: The theory possess l^"/^ x SU{nf), where the first factor is the enhanced 
topological times discrete xS. A: All symmetries broken. B: l}/^ broken. C: None of the symmetries 
is broken. D: Discrete and continuous xS broken. In this work, we study the transition from B to C 
by reliable semi-classical field-theoretic methods. The thick green line in the lower left corner is the 
phase boundary at small L where the theory is under complete analytic control. The field theory in 
the vicinity of this fine can microscopically be mapped, for SU(2) and S'?7(2)/Z2, respectively, to two 
types of XY-spin systems with perturbations. The dashed lines represent the conjectured extrapolation 
where the different phases coexist. 

The phase diagram of the SU{2) theory in the L — 13 plane (here, /3 = ^) is shown in 
Fig. 1. According to our current understanding of QCD(adj), the theory should have four 
phases. Phase-A exhibits confinement with discrete and continuous chiral symmetry breaking 
(xSB). Phase-B has confinement with discrete but without continuous xSB. The existence 
of this phase can be shown analytically [8]. Phase-C is a deconfined, chirally symmetric 
phase. The existence of this phase can also be shown analytically [23] . Phase-D is deconfined 
with discrete and continuous xSB. The existence of this phase is understood numerically 
through the lattice studies [24]. Our current work addresses, by reliable continuum field 
theory techniques, the second-order transition between B and C. 

A few comments on the B-C phase boundary are now due. In different theories, this 
phase boundary has different interpretation. 

1. In the SU{2)/'L2 theory as well as in the Z4-spin system, the Z4 is broken in the low 
temperature phase and restored at high temperature phase. 

2. In the SU{2) theory as well as in the associated spin system, the symmetry is Z2 x Zg"^. 
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It is broken down to 7L\ in the low-temperature phase, i.e., center symmetry is unbroken. 
At high temperatures, the symmetry is broken down to Zg^. 

This difference is there because the set of electric and magnetic charges that we are allowed 
to probe the two gauge theories and spin systems are different. 

1.2 The theory of melting of 2d-crystals and SU(3) QCD(adj) 

Our next result concerns ^{/(Sj/Za QCD(adj) with n/ massless adjoint Weyl fermions on 
M^'^ X S];^. The theory near the deconfinement transition is described by a spin model, which 
is a "vector" generalization of (1.1). This is the theory of two coupled XY-spins, described 
by two compact variables whose periodicity is determined by the S\J (3) root lattice: 

= {el, el) = e^ + 2^1 = 4 + 2^2 . (1.5) 

Here di {i = 1,2) are the simple roots of SU{3), which can be taken to be ai=(l,0), 
a2=(— 5, ■^). The theory is defined, similar to (1.1), by a lattice partition function with: 

3 _ 3 

x;fi=l,2 i=l X i=l 

where Vi {i = 1,2,3) are the weights of the defining representation (see Section 4.3 for their 
explicit definition) and the sum in the potential term includes also the affine root. 

The first term in (1.6) is essentially^ the model used in [5] to describe the melting of 
a two-dimensional crystal with a triangular lattice. There, the two-vector 0^. = (^^,0^) 
parametrizes fluctuations of the positions of atoms in a 2d crystal around equilibrium (the 
distortion field) and n is proportional to the Lame coefficients. A vortex of the compact vector 
field e^ describes a dislocation in the 2d crystal. The Burgers' vector of the dislocation is the 
winding number of the vortex, now also a two-component vector. Without going to details 
of the order parameters and phase diagram, see [5], we only mention that the melting of the 
crystal occurs due to the proliferation of dislocations at high temperature, which destroys the 
algebraic long-range translational order. 

To relate (1.6) to our theory of interest, SU{2>) QCD(adj), we shall argue in Section 
4.3, that the vortices in (1.6) describe the electric excitations, the VF-bosons, in the thermal 
theory. The vortex-vortex coupling k is related to the parameters of the underlying theory 
exactly as in (1.2). In the melting applications of the "vector" XY-model, there is an exact 
C/(l) X U{\) global symmetry e^^ex + c. This symmetry forbids terms in the action which 
are not periodic functions of the difference operator. In QCD(adj), it is broken to Z3 x Z3 by 
the potential term in (1.6). This breaking is crucial and in the chosen duality frame describes 
the effect of the magnetically charged particles in the thermal gauge theory — the magnetic 
bions responsible for confinement. The two Z3 symmetries of (1.6) are the topological Z3 
(associated with the nontrivial 7ri(S'?7(3)/Z3)) and the Z3 discrete subgroup of the chiral 
symmetry of the theory. 

^In disguise: see eqn. (4.30) for another description of (1.6), making its relation to [5] more obvious. 
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The nature of the deconfinement transition in the SU{3) QCD(adj) is not yet completely 
understood. One obvious feature that follows from (1.6) is that at low temperature (large 
K, see (1.2)), where vortices can be neglected, the topological/chiral Z3 x Z3 symmetry is 
spontaneously broken, as expected in the confining phase (and also follows from the zero- 
temperature analysis). 

The SU{3) theory also exhibits electric magnetic duality, which interchanges electric and 
magnetic fugacities and inverts the coupling in (1.6): 

(1.7) 

K 

2 

Thus, a candidate for the deconfinement transition is the self-dual point Tc = ^ . At 
present, we do not know if at the self-dual point (1.6) is a CFT (which would be the case if 
the transition was continuous). 

In Section 4.2, we show that the leading order (in fugacities) renormalization group 
equation for k does, indeed, exhibit a fixed point at the self-dual point of e-m duality (1.7). 
However, both 1^-boson and bion fugacities are relevant at that point — indicating that a 
strong-coupling analysis is necessary, even if the fugacities are small at the lattice cutoff, 
unlike the SU{2) case. We note that in SU{3) QCD(adj) it should also be possible to use 
CFT methods at the self-dual point, as in [4], to study the existence and nature of the critical 
theory. 

1.3 Remarks on the afRne XY-spin model and SU(Nc > 3) 

We define the affine XY-spin model by: 

N 

(_^^)affine ^ ^ S ^l?. • (4+A " 0.) , 9. = C + ^TTOi . (1.8) 

x;fi=l,2 i=l 

This is a natural Lie-algebraic generalization of the XY model, and possesses a U{1)^''~^ 
global symmetry, 6^ 6x + c, where 6 is now an (A'^c—l) -dimensional vector with periodicity 
defined by the simple roots of SU{Nc)- The physics of the phase transition in this model is 
expected to be a generalization of BKT. However, as in the discussion of SU{2) and SU{3) 
theories, the spin systems relevant for gauge theories contain certain symmetry-breaking 
perturbations. There are two such interesting deformations of (1.8) which relate it to gauge 
theories. One preserves a Zjv subgroup of the global symmetry and the other preserves a 
Zjv X Z^r symmetry (which accidentally enhances to Z4 for N = 2). These deformations are:^ 



-I3H = {-PH)''^''^ + { 



Ex Ej=i y cos 2ai ■ 6^ 



7t 



Ex Eili y cos 2 (a, - Qi_i) • 4 , X z*^^ 



(1.9) 

N ' 



^As noted earlier, there is a well-known classification of phases for the (7(1) — >■ Zp model, which crucially 
depend on p, see e.g., [22]. It is crucial to note that in U(l)^~^ — > Zn, N does not play the same role as p, 
and the classification is different. 
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where the sum now includes also the affine root (see Section 4.1). 

As in the discussion of SU (2) and SU{3), the vortices of the affine model can be identified 
as T4^-bosons in the gauge theory, forming an electric plasma with interactions dictated by 
the Cartan matrix • Oj. The perturbations correspond to the proliferation of charges, 
corresponding to magnetic monopoles for the case and magnetic bions for x Z^ case. 
The interactions in the case of magnetic monopoles are also governed by the Cartan matrix, 
up to an overall coupling inversion.^ 

The affine Zjv-model is appropriate for a bosonic gauge theory, such as deformed Yang- 
Mills theory [25], or equivalently, massive QCD(adj) in a regime of sufficiently small L. 
The Zjv symmetry in the spin-system is associated with the topological Z^ in the bosonic 
SU{N) /'Ln theory. The vortex-charge system exhibits e-m duality. The charges (magnetic 
monopoles) and vortices (VF-bosons) in the plasma have the same interactions, dictated by 
the Cartan matrix Oj • dj for both e-e and m-m interactions. 

There are only a few remarks about iVc > 3 QCD(adj) that we will make in this paper. 
First, the Ztv x Zjv symmetry in (1.9) is associated with the topological Z^ and discrete Z^ 
chiral symmetries of SU{N)/Xiy QCD(adj). We show in Section 4.1 that the SU{Nc > 3) 
QCD(adj) theory on M^'^ x Sj^ near the deconfinement temperature is described by an e-m 
Coulomb gas which does not exhibit e-m duality. This is in contrast to the Nc > 3 case of 
the deformed Yang-Mills theory, or massive QCD(adj). The absence of duality in the case of 
QCD(adj) is due to the fact that, for Nc > 3, the magnetic (bions) and electric (VF-bosons) 
particles in the gas have different interactions, governed by: 

di ■ dj = 6ij - - for W - bosons , 

Qi ■ Qj = 36ij - 25ij+i - 25jj_i + + fo'^ magnetic bions . (1.12) 



*A brief note on the symmetries of (1.9) is due. The Zjv symmetry acts on 9 by shifts on the weight lattice, 
Zjv : 9 — > 9 — 2'Kkvi, k = 0, . . . N ~ 1, for any i, 

g2»Pj-'' ^ g-47rfeiiA,-.Pig2i,7j.e _ ^i2^^2ivj-e ^ (1.10) 

where we used identities given in Footnote 29 in Section 4.3 (it might appear that there are — 1 symmetries 
in (1.10), but note that the difference between transformations with different i's is a shift of by a root 
vector, which by (1.8) is an identification, not a symmetry). Since the "monopole" perturbation in (1.9) 
is e^'"j * = e^"^J *e~^"^J+i", its invariance under Z]^ is manifest. The Z^v invariance then trivially holds 
also for bion perturbations g^^aj Sg-aiaj^i-e^ rpj^^ bion-induced perturbation (the one on the bottom line in 
(1.9)) has an extra independent symmetry, which rotates the monopole operator by N^^ root of unity, 
^2iaj-9 _^ g22c«j -9 ^^i^^ Symmetry is a remnant of the discrete chiral symmetry in the gauge theory. In 
order to see the action of the symmetry on the 9 field, we define the Weyl vector, p = X]i=i Mm where jli 
axe the fundamental weights of SU{N), defined through the reciprocity relation di ■ fij = Note that p 

satisfies p- di = ^ for i = 1, . . . , N — 1, and p ■ dn = Consequently, 



Z^^: S^S+'^-^p, q = 0,...N 

. 47rg - 



iV 

3j-p^2iaj-tt _ ^i^^zicj-a ^ (1.11) 

Therefore, the bion induced term in (1.9) is also invariant under the Z^ symmetry. 
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The magnetic bions (of charges Qi = di — Oi-i) have next-to- nearest-neighbor interactions, in 
the sense of their "location" on the extended Dynkin diagram, while the Vl^-bosons only have 
nearest-neighbor interactions. The different interactions are ultimately due to the particular 
composite nature of bions, which is itself is due to the presence of massless fermions in the 
theory. 

In Section 4.2 and Appendix C, we derive the renormalization group equations for the 
SU{Nc) QCD(adj) Coulomb gas to leading order in the fugacities and show that none of them 
displays a fixed point. Thus, the picture of the deconfinement transition for > 3 is, so far, 
not complete. 

In Section 5, we list possible avenues for further studies. 

2. The zero temperature dynamics of QCD(adj) on M}'"^ x 
2.1 Review of perturbative dynamics at zero temperature 

We consider four-dimensional (4d) SU (Nc) Yang-Mills theory with nj massless Weyl fermions 
in the adjoint representation, a class of QCD-like (vector) theories usually denoted QCD(adj). 
The action for SU{Nc) QCD(adj) defined on M^'^ x is: 



S= tr 



1 



2g^ 



(2.1) 



where we use the signature r]^^, = diag(+, — , — , — ), Fmn is the field strength, D]\f is the 
covariant derivative, / is the flavor index. A/ = Xi^ata, o, = 1,...,N^ — 1, ctm = (1,'r), 
o'M = (1)~^)) 1 and r are respectively the identity and Pauli matrices, and the hermitean 
generators ta, taken in the fundamental representation, are normalized as Tvtatb = 
The upper case Latin letters M,N run over 0,1,2,3, while the Greek letters /i, run over 
0, 1, 2. The components and 3 denote time and the compact dimension, respectively. Thus 
= + L, where L is the circumference of the circle. 
Notice that we are compactifying a spatial direction and the time direction £ M^'^ is non- 
compact. Thus the fermions and the gauge fields obey periodic boundary conditions around 
S]^. The action (2.1) has a classical U{1) x SU{nf) chiral symmetry acting on the fermions. 
The quantum theory has a dynamical strong scale Aqcd such that, to one-loop order, we 
have: 

9/ , IGtt^ 1 

^ ^ = ^ 1 ( 2/A2 ^ ' 2.2 

/3o log(/i7A^^^) 

where fi is the renormalization scale and /3o = ^^^^ ^^i^" ; we will only consider a small 
number of Weyl adjoint fermions, nj < 5.5, to preserve asymptotic freedom. The theory 
also has BPST instanton solutions with 2A'^c zero modes for every Weyl fermion and the 

N 

corresponding 't Hooft vertex, ~ e (deti^jA/^aA j^a) explicitly breaks the classical U{1) 
chiral symmetry to its anomaly-free '^2Ncnf subgroup. The classical SU{nf) x U{1) chiral 
symmetry of the theory is thus reduced to: 



SU{nf) X 'L2N^nf 



(2.3) 
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where the factor common to SU{nf) and 1j2Ncnf is factored out to prevent double count- 
ing. The Z2 subgroup of the 1j2Ncnf is fermion number modulo 2, (—1)^, which cannot be 
spontaneously broken so long as Lorentz symmetry is unbroken. Thus, the only genuine dis- 
crete chiral symmetry of SU{Nc) QCD(adj) which may potentially be broken is the remaining 
Zjv^, which we label as Zi^, irrespective of the number of flavors. 

An important variable is the Polyakov loop, or holonomy, defined as the path ordered 
exponent in the direction: 

nL{x) = Pe-^'i^'^"'"'\ (2.4) 

where x £ M^'^. This quantity transforms under x-dependent gauge transformation as — )• 
U^^^}lU, and hence the set of eigenvalues of fi^, is gauge invariant. The gauge invariant trace 
of the Polyakov loop trQ^ serves as an order parameter for the global center symmetry 

j 2TTk 

trr^i — 7- e trO,L, with k = 1, . . . Nc- If we take LAqqd ^ 1, a reliable one-loop analysis 
can be performed to study the realization of the zj^^^ center symmetry. Integrating out the 
heavy Kaluza-Klein modes along S^, with periodic boundary conditions for the fermions — 
remembering that our is a spatial, not a thermal circle — results in an effective potential 
for Q: 

2 °° 1 

Vesm = (-1 + ^/)^ E -J^'^ii'' ■ (2-5) 

n=l 

Since 0^,(2;) = e^^'^ai^) ^ ggg gqj^_ (2.4), the action for the x'^-independent modes of the gauge 
field is:^ 



Sn = I ^tr 

L,2 



-F^uF''- + {D^A^f - ^V,ff{As) + 2ig'y {a'^D,, - ta^ [A3,]) A, 



(2.6) 



The minimum of the potential V^fj (2.5) for ?t,j > 1 is located, up to conjugation by gauge 
transformations, at [7]: 

(J^l) = diag(l,WAr^,w^^, . . .wj^;;"^) , ujN, = e^- , (2.7) 

where ry = e^c for even Nc and 77 = 1 otherwise. Since tr(0) = 0, the Z]^^ center symmetry 
is preserved. This is in contrast with the nj = case, where the theory on §^ is equivalent to 
a finite temperature pure Yang-Mills theory, where the theory deconfines at sufficiently high 
temperature [23]. In the supersymmetric nj = 1 case the one-loop potential vanishes and the 
center-symmetric vacuum is stabilized due to non-perturbative corrections [26-28]. 

In the vacuum (2.7), the SU{Nc) gauge symmetry is broken^ by the Higgs mechanism 
down to C/(l)^'^^^. Because the unitary "Higgs field" JIl is in the adjoint, its off-diagonal 

""To avoid confusion, we note that in subsequent sections we study the Euclidean theory, where, for conve- 
nience, we relabel the compact direction x'*' and the "Higgs field" A4,. 

®The dynamical abelianization on M}''^ x S\ in QCD(adj), as in the Coulomb branch of the Seiberg-Witten 
theory [26] and the 3d Polyakov model [29], is the ultimate reason that the theory admits a semi-classical 
analysis. It also helps to continuously connect two confinement mechanisms, which were previously thought 
to be unrelated [30]. See [31] for a bosonic version. Another interesting recent work taking advantage of the 
abelianization, this time on R^'^ x S^, is [32]. 
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components are eaten by the gauge fields which become massive. At the lowest mass level, the 
T^-bosons have mass mw = yn~', the heavier VK-bosons have mass with k = 2,. . . Nc- 
The Nc — 1 diagonal components of the Higgs field also obtain masses ~ y^g^Nc, due 
to the one-loop potential ^(^3) [33]. In addition, the off-diagonal fermion components (the 
color space components of A that do not commute with {A3) also acquire mass, with masses 
identical to those of the VF-bosons). The coupling g ceases to run at 1/L ^ ANc- 

The conclusion is, then, that for distances greater than LNc (not L), the three-dimensional 
(3d) low energy theory is described in terms of the diagonal components of the gauge field 
(the Nc — 1 "photons") and the {Nc — 1) x Uf diagonal components of the adjoint fermions 
(denoted by Xk,i, with 1 = 1,... nj, k = 1, . . . Nc — 1). At the perturbative level, up to higher 
dimensional operators suppressed by powers of L ~ , the Lagrangian is: 

^ ^ Nc-l N-1 «/ 

S= -Ao^T. Ph.^ + E E ' (2-8) 

■^its ' y ■ — 1 7 — 1 r 1 



k=l 1=1 



i.e., that of a free field theory. In the next section, we review how non-perturbative effects 
change this Lagrangian. 

2.2 Review of the SU(2) non-perturbative dynamics at zero temperature 

We will specialize to A'c = 2 in this section. To further cast the free theory (2.8) into a 
form that will be useful in our further analysis, we note that we can dualize the 3d "photon" 
into a periodic scalar field — the "dual photon." The duality transformation can be performed 
by first introducing a Lagrange multiplier field a imposing the Bianchi identity. Thus, we 
introduce an additional term to the action S, i.e., eqn. (2.8) with Nc = 2: 



6S=^ [ ae^uxd'^F'^^ 



(2.9) 



Note that we are still considering the duality in Minkowskian signature and that in Euclidean 
space, a factor of i will appear. Integrating over F^i, in the path integral then amounts to 
substituting: 

F^, = -^fM-A^V (2.10) 
into S + 6S. The result is the action written in terms of the "dual photon" field a: 



1=1 

which captures the long-distance physics of the photons to all orders in perturbation theory 
and to leading order in the derivative expansion. 

The next step is to consider the modification of (2.11) by non-perturbative effects. How- 
ever, before we study the non-perturbative dynamics, we shall make some comments on SU (2) 
vs. SU{2)/7j2 theory, discrete fluxes, and the topological Z2 symmetry. Our comments are 
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complementary to those of [34]; they also generahze to higher-rank groups and wih be pre- 
sented elsewhere in more detail. The brief discussion in Section 2.2.1 below should be useful 
to get a picture of the origin of the topological Z2 symmetry — which plays an important role 
in the analysis of the finite-temperature theory and the deconfinement transition in SU (2) /Z2 
theory. For SU{2), the center symmetry plays an analogous role. 

2.2.1 Hamiltonian quantization of compact U(l), SO(3) vs. SU(2), discrete fluxes, 
and the topological Z2 symmetry. 

In this section, we explain the difference between the SU{2) vs. 5C/(2)/Z2 theory in the 
Hamiltonian formulation. Let us first focus on the bosonic part of (2.6), without considering 
fermions. Such a bosonic theory arises in the study of QCD(adj) with massive fermions or 
deformed Yang-Mills theory. This bosonic theory, in the long-distance regime and to all-orders 
in perturbation theory, admits a description in terms of the bosonic part of (2.11): 

JmixT2 2-L V47ry 

where we are now considering the theory with space additionally compactified on a large 
of area Aj2, i.e., on x T^. This is helpful to understand the appearance of discrete fluxes 
and the associated topological Z2 symmetry. 

It is well known that on T^, a compact U{1) theory has magnetic flux sectors, see, 
e.g., [22], labeled by integers n counting the unit of flux quanta on and that, for a large 
T2, the energies of these flux sectors are almost degenerate. To see these sectors in the dual- 
photon language, consider the dynamics of the zero mode ctq of the a- field on T^. The zero 
mode is governed by the quantum-mechanical action: 

The action 5*0 is obtained upon integration of (2.12) over the by keeping only the zero 
mode. Clearly, the zero mode dynamics is that of a free "particle" of "mass" M. To see that 
the "particle" is actually moving on a circle and infer the size of the circle, i.e., its radius R, 
o'o = f^o + ^ttR, consider the solution of the equation of motion, ao = t V + const, where V 
is the constant "velocity" of the particle. From the duality relation (2.10), a constant value 
oi ao = V implies that there is a constant magnetic field on the T^: 

F12 = -/y V . (2.14) 

At this point, it is useful to ask whether the nonabelian theory underlying the U{1) is an 
S'0(3)=S'C/(2)/Z2 or an SU{2) theory. These are distinguished by the kinds of C/(l)-charges 
that are allowed in the spectrum, as we now explain: 

1. In the 5'0(3) theory, a charged particle can only have integer charge q = ±1, ... under 
the unbroken U{1). In order that the phase factor e^'^^c^as^ = e^'^^s^ acquired by 
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such a particle upon its motion on be unambiguously defined,^ it must be that 
Jj2 B = Fi2Aj2 = 27rn, with n G Z. Eqn. (2.14) then implies that the "velocity" of ctq 
is quantized, V = n ^^l"" '^^^ momentum conjugate to do is IIo-o = MV = ^ and the 
energy is: 



2M 2M \2J ' [n^oEZ, for n G 22 

We refer to the sector with half-integer (integer) conjugate momenta IIo-o as half-integer 
(integer) flux sector. The expression of the energy (2.15) is the same obtained upon 
quantizing the motion of a particle of mass M on a circle of radius R = 2. This implies 
that the coordinate do of the particle from (2.13) is a periodic variable of period 47r. 
Hence, we conclude that the dual photon a in the SO{'i) theory is a compact field with 
periodicity 2tt: 

o- = cr + 47r , for the 50(3) theory. (2.16) 

Furthermore, since M ~ ^^2, see (2.13), the above-mentioned near-degeneracy of the 
flux sectors at large is evident from (2.15). 

As functions of the zero mode (Tq of the dual photon on T^, the flux sectors have 
Schrodinger wave functionals ^'„,((Jo) = e*"'^"/^. We shall later argue that shifts of a 
by 271 play a special role, namely that such shifts constitute a symmetry of the «S'0(3) 
theory, called the "topological Z2 symmetry," associated with the nontrivial 7ri(50(3)). 
Under such shifts of cr, we have: 

Z*: do ^ do + 2^, ^„(ao) ^ (-)"^'„(do). (2.17) 

All states in the Hilbert space (at finite T^) of the 50(3) theory are either even or odd 
under the topological Z2, as in (2.17). 

2. In the SU{2) theory, on the other hand, particles can have half-integer charge q = ib^, ... 
under the U{1) (imagine adding a heavy "electron" in the fundamental representation to 
QCD(adj)). In order that their propagation on be well defined (i.e., the phase factor 
6*2/0"^ be unambiguous), it must be that ^Fi2Aj2 = 27m, n G Z. Then, the "velocity" 
of do is quantized as V = 2n ^ , the conjugate momentum is Ilo-g = MV = n, and 
the energy is: 

^'""^'^ = = Jm""' ' "'^^ ^' for n e Z . (2.18) 

Thus, in the SU{2) theory, only the integer flux sector exists. Eqn. (2.18) is exactly 
the quantum of energy of a particle of mass M moving on a circle of radius R = 1. 

^Consider an oriented Wilson loop C on a a compact two dimensional space such as T^. There is no 
unique surface whose boundary is C. Both "inside" and "outside" surfaces, Si and So, are acceptable choices 
leading to a two-fold ambiguity, which is resolved provided e'^-'^c?'* = e'^^^^ ^ _ Js^ b ^ j ^-^^ phase is 
independent of the choice. This implies, q B + q B = q J^2 B — 2Tvn. 
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Thus, the coordinate of the particle do in the SU{2) theory has period 2-kR = 27r. We 
conclude that the dual photon in the SU{2) theory is a compact field, but with half the 
periodicity of (2.16): 

fj = cr + 27r , for the SU{2) theory. (2.19) 

Next, we recall that the compact U{\) theory (2.12) is a low-energy remnant of the broken 
nonabelian theory and that the latter has monopole-instanton solutions (at sufficiently large 
one can use the infinite volume solutions). The monopole-instantons are finite Euclidean 
action solutions and play an important role in the non- perturb at ive dynamics. They affect 
the perturbative long-distance Lagrangian (2.11) in a manner to be described in the following 
section. 

These solutions are constructed from the 4d static 't Hooft-Polyakov monopoles as re- 
cently reviewed, e.g., in [35]. As explained there, in a U{1) theory which descends from an 
SU(2) or 50(3) theory with a compact adjoint Higgs field, there are two types of monopoles,^ 
M.I and M.2-, and their anti-monopoles, see also Section 2.2.2 (the second type of monopole 
is present due to the topology of the adjoint Higgs-field). However, most of the argument 
in the bosonic theory described below is independent of this aspect, because the magnetic 
charge and action of A^i and A^2 (not M2) coincide. These classical solutions are, of course, 
identical in the SU{2) and 5*0(3) theories. The magnetic charge of a 't Hooft-Polyakov 
monopole-instanton is: 

/ B = \ f e^.xF"^ = inQm, with \Qm\ = 1 , (2.20) 

i.e., we take the abelian long-range field of a monopole to be i? = Qm^, without a factor of 
4tt in the denominator. The Dirac quantization condition is [45]: 

QmQe = ^, neZ, (2.21) 

and the electric charge of VF-bosons is Qe = ±1. The 't Hooft-Polyakov monopole solutions 
have Qm = =tl- This is twice the allowed minimal value of magnetic charge in the 50(3) 
theory with Q™" = il, and is exactly equal to the minimal value of magnetic charge in the 
SU{2) theory, where Q™" = ±1 instead. 

The Hilbert space interpretation of the monopole-instantons is that they facilitate tunnel- 
ing between states of different magnetic ffux on T^. From (2.20), it is clear that a monopole- 
instanton tunneling event changes the magnetic flux on by one unit of magnetic flux.^ 
We note, for completeness, that the tunneling events that connect the vacua in different 



Of lowest action, which are most relevant in the small-L regime we study here; there is an infinite Kaluza- 
Klein tower of monopole-instanton solutions of increasing action. 

^To see this, deform the integral (2.20) to an integral over two infinite planes separated by a large 
Euclidean time interval 2T, i.e., An = J^2 B = J^2 t It^ ^''^'^ recall that a "half-unit" flux on 

corresponds to B = 2tv; to avoid confusion, recall that after eqn. (2.15) we adopted the convention of 
calling the minimal flux allowed "half-unit", corresponding to a half-integer n in eqn. (2.22) below. 
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SO(3) Flux-states and tunneling events 
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SU(2) Flux-states and tunneling events 
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Figure 2: In an SO{3) theory on T^, the flux sectors are associated with integers and half-integers, 
as per eqn. (2.15). In an SU{2) theory, only integer flux sectors exist, as (2.18) shows. In both, the 
tunneling events (which are a realization of the 't Hooft-Polyakov monopole-instanton on T^) change 
flux by one unit. Consequently, the SU{2) theory has a unique vacuum, which is a superposition of 
all flux sectors, whereas the 5*0(3) theory has two vacua, distinguished by their Z2 quantum number 
(2.17): one is the superposition of half-integer flux sectors and the other is the one of the integer flux 
sectors. This is the Hamiltonian realization of topological Z2 symmetry. 

flux sectors are of two- types, call them Ii and I2, these are the descendants of TWi and M2 
monopoles mentioned above, and they lead to change in the flux by one unit, i.e., 



For the purpose of our discussion, Ii and I2 play the same role as both increase the flux by 
one unit.^^ To study the consequences of the tunneling between the various flux sectors, we 
have to, again, make a distinction between the S0{3) and SU{2) theories (as illustrated on 
Figure 2): 

1. In the SU{2) theory, the 't Hooft-Polyakov instanton-monopoles facilitate tunneling 
between sectors of magnetic flux differing by one unit of flux quantum. Since the SU (2) 
theory only allows an integer number of magnetic flux quanta, as per (2.18), one expects 
that at sufficiently large the true ground state is a unique superposition of all flux 
sectors. 

2. In the SO{2>) theory, both integer and half-integer sectors are allowed, (2.15), but, as 
in the SU{2) theory, tunneling events can only change the flux on by one unit of 

^°In the Hamiltonian formulation of the Polyakov model, we would only have the Ji-type instanton. This 
follows from the non-compactness of the adjoint Higgs. In deformed Yang-Mills, or massive QCD(adj), due 
to the compact nature of the adjoint Higgs, there are two-types of flux changing instantons. Also, note that 
the existence of this second type of instanton is the main difference with respect to ref. [22]. The existence 
of the two-types of flux-changing instantons has important physical implications in (at least) two cases: when 
massless fermions or a non vanishing 6 angle is introduced. These aspects will be discussed elsewhere. 






(2.22) 
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flux — there is only tunneling between sectors of integer flux and between sectors of half- 
integer flux, but not between integer and half-integer flux sectors. Thus, the broken 
SO{3) theory actually has two flux sectors that do not interact with each other. These 
can be labeled by a Z2 quantum number, as in (2.17): states in Hilbert space have 
eigenvalue +1 under the topological Z2 symmetry if they are in the integer flux sector 
and —1 if they are in the half-integer flux sector. 

The integer and half-integer flux sectors in the U{1) theory are the long-distance rem- 
nants of the two 't Hooft magnetic flux sectors^'^ [36] in SO{3) theories on T^. In the 
limit of large they become degenerate. The even and odd magnetic flux^^ sectors 
in 50(3) theories on are distinct superselection sectors of the theory, unless local 
operators on the Hilbert space that can change the value of the Z2 flux are introduced. 

The topological Z2 symmetry (2.17) is generated by an infinitely large spacelike Wilson 
loop [38]. Let C denote boundary of a surface S, which we eventually take to be M^, 
formally, C = dM?. In the low-energy theory of the dual photon, using (2.10) and the 
definition of conjugate momenta, the Wilson loop operator can be written as: 

W = e^^c^'^'- = ^ifBtfix ^ ^i2-K J (fxiia{x) ^ (2.23) 

where Ii(j{x) is the momentum canonically conjugate to a, i.e., with equal time commu- 
tator [no-(x), (T(y)] = —iS{x — y) (the equal time argument is suppressed). The physical 
interpretation of W is that it measures the amount of Z2 magnetic fiux present at the 
given time. The operator e*°"(^)/^ and W obey the equal time commutation relation: 

^gia(x)/2 ^ _gia{x)/2^^ ^2.24) 

which is the 3d version of 't Hooft algebra for 50(3). The interpretation of the oper- 
ator 6*°"^^)/^ is that it creates a unit of Z2 flux at x (a Z2 vortex). In the context of 
the underlying M^'^ x S^]^ theory, the commutation relation (2.23) may be seen as the 
dimensional reduction of the 't Hooft commutation relation: 

W{C)T{C') = e*"^(^'^')r(C')iy(C7) (2.25) 

where i{C, C) is the linking number of the two closed loops C and C , and T{C') is the 
't Hooft operator. We consider the limit in which the inside of C fills M?. e*°^^^'-*/^ is the 
long-distance "residue" of a 't Hooft loop of a Q™" = | monopole in the 4d theory, 
winding around and located at x E M^, hence, i{C, C) = 1. 

If we allow the "Z2-vortex" operator e'^'^i^)/'^ in the Hilbert space, see [33], correlation 
functions like (e»'^(^)/2e-*'^(y)/2) 

are observable, and we have to consider the possibility 

^'^In the full theory, the Z2-odd ("half-integer" in our convention) flux sector is constructed by considering 
SO{'i) bundles on twisted along one of the non contractible loops on by gauge transformations in the 
topologically nontrivial class 7ri(S'0(3)) = Z2. This way to construct discrete flux sectors on tori is explained 
in [36]; see [37] for different perspectives. 

^^Also their T'^ generalization and the related discrete electric flux sectors, which we do not consider here. 
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of spontaneous symmetry breaking of the topological Z2 symmetry. Thus, in the limit 
of large T^, the ground state of the Hamiltonian may not be an eigenstate of Z2 as in 
(2.17), but may be a linear superposition between even and odd sectors corresponding 
to a fixed value of (e*°"/^) = ±1 instead; we shall see in the next section that this is 
indeed the case in the SO {3) theory. 

The remarks in this section concerning the topological Z2 symmetry and the difference be- 
tween bosonic SU(2) and 50(3) theories will be useful when considering the different de- 
scriptions of the deconfinement transition in Section 3.3.3. 

Fermionic zero modes and flux sectors: The above discussion refers to a theory without 
fermions. In the theory with fermions, because of fermion zero modes, tunneling events that 
change the magnetic flux |n) — t- |n + l) (n = 0,ibl/2,ibl) cannot occur without also changing 
the fermion number of the state. However, tunneling events between states of zero fermion 
number that change the flux sectors as \n) — t- |n + 2), do occur. This effect is a descendant 
of magnetic bions, a certain type of topological molecule which changes the magnetic flux by 
two units (their magnetic charge is two times the 't Hooft-Polyakov monopole charge or four 
times the | -vortex charge mentioned above) and that have no fermion zero modes, see Figure 
3. These defects on M'^ x S-*^ will be reviewed in the next Section 2.2.2. Thus, in the SU{2) 
theory with fermions we will obtain two ground states, almost degenerate at sufficiently large 
2d volume, while in the SO (3) theory, there will be four almost degenerate ground states. 
The detailed structure of the ground states in theories with fermions can be given via the 
quantum mechanics of zero modes of both the dual photon and the massless components of 
the fermions and will be discussed elsewhere. 

To sum up this section, in the bosonic theory, SU{2) has a unique vacuum, and trivial 
topological symmetry, SU (2) /7,2 has two vacua and a Z2 topological symmetry. The pres- 
ence of adjoint fermions doubles the number of vacua in each case, making it two and four, 
respectively. 

2.2.2 Nonperturbative eff'ects in the SU(2) theory on M^'^ x §1. 

The perturbative Lagrangian (2.8) in terms of the compact dual photons and neutral fermions 
does not properly reflect the dynamics of the long-distance theory. There are important non- 
perturbative effects due to the existence of monopole-instantons in the broken QCD(adj) 
theory. These are not only qualitatively, but also quantitatively important, and their con- 
tributions are under full theoretical control in the small-L weak-coupling regime. In what 
follows, we will briefly describe the properties of the relevant monopole-instantons and their 
role in the dynamics. 

In this section, we will, for the most part, have in mind the dual photon appropriate to 
the SU{2) theory, i.e., a with periodicity (2.19). We stress that all previous studies of gauge 
theories on M"^ x (supersymmetric or otherwise) have been also in the context of SU (2) 
theories. 

^^If long-distance correlations in the values of e^"^^ are present, not allowing for Z2 symmetry-breaking 
states leads to violation of cluster decomposition, see, e.g., Ch. 23 in [39]. 
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^0 = o X = -X- (2-26) 



Because the effective 3d theory (2.11) originates in a 4d theory compactified on M'^ x S]^ 
(we now consider the Euchdean setting appropriate for studying instantons), there is a whole 
Kaluza-Klein tower of monopole-instantons, in addition to the single fundamental monopole- 
instanton that exists in a purely 3d theory. These solutions have been extensively described 
in the literature [40,41]. In this paper, we work in the small-L regime and to leading order 
in the semi-classical expansion — so that the contributions of the higher-action Kaluza-Klein 
monopole-instantons will not be important (unlike the study of ref. [30]). Thus, we will only 
use the lowest action monopole-instantons. Below, we will describe the two types of solutions 
of lowest action of the self-duality equations Fmn = Fmn = \(-mnklF^^ ■ In the center- 
symmetric vacuum (2.7), their action equals \ the BPST instanton action. We denote the 
action by Sq: 

1 87r2 _ 47r2 

2 ""^"^ 

Because the solutions are self-dual, their topological charges, Qt = /r3x§i ^mn^mn^ 
are all equal to each other and to | in the center-symmetric vacuum. The unbroken-f/(l) 
magnetic field far from the core of the two types of self-dual monopole-instantons is given by 
(for a review of the explicit solutions, see, e.g., [35]): 

= ±^ . (2.27) 

Thus, according to our normalization (2.20) the two lowest-action solutions — the monopole- 
instanton (+) and the KK-monopole-instanton (— ) (as the extra self-dual solution is often 
called) — have magnetic charge ±1. We stress again that the + and — charge monopole- 
instanton solutions are both self-dual and that the corresponding anti-self-dual antimonopole 
solutions carry opposite magnetic charges. 

In a theory with massless adjoint fermions, the monopole-instantons have fermionic zero 
modes (see [42,43] for the relevant index theorem) and generate, instead of the usual monopole 
operator e~^°e^^'^ , operators of the form: 

Ml = e-^oe^'^ dej; AjA/ , M2 = e'^^e-'" det A jA/ , (2.28) 

where the determinant is over the n/ flavor indices. The form of the monopole-instanton 
induced operators (2.28) has interesting implications for the physics of mass gap generation 
and confinement: 

• First, sincG und-Gr the "^Anf 

discrete chiral symmetry of (2.3) detjj XjXi — )• — det/_j AjA/, 
the invariance of the monopole-instanton amplitude under the exact discrete chiral sym- 
metry demands that the dual photon transform as well, as shown on the second line 
below: 



Z'^'^ : det AjA/ ^ - det AjA/ 



(2.29) 



- 18 - 



Here, TL,^ denotes the genuine discrete chiral symmetry which cannot be rotated away 
by a discrete nonabehan flavor transformation, or be included in (—1)^, as explained 
just after eqn. (2.3). Thus, e^" is the topological order parameter (disorder operator) 
associated with the Zg'^ symmetry o" — )• a + vr. 

• Second, symmetry considerations alone show that the 1^ symmetry permits a purely 
bosonic potential term, ~cos2(7, in the long-distance effective action. This term was 
shown [8, 9] to be due to a novel type of topological molecule and referred to as "mag- 
netic bions." In a Euclidean context, where monopoles are viewed as classical particles, 
the magnetic bions are molecular (correlated) instanton events of self-dual monopole- 
instantons and anti-self-dual KK-monopole-instantons that carry twice the charge of a 
magnetic monopole. Due to this reason, we sometimes refer to these defects as topologi- 
cal molecules. Magnetic bions arise in second order in the semi-classical e"'^" expansion 
as ~ e~^'^° cos 2o". 



Considering the second point in some more detail, we note that the "magnetic bion" topolog- 
ical molecules are stable "bound states," where the magnetic repulsion between the charge 
monopole and the charge anti-KK- monopole is balanced by attraction due to fermion 
exchange, resulting in an interaction potential of the form: 



^bion(?') 



h 4n f In r , 



(2.30) 



with stabilization radius r* = -KL/g^nf. The resulting "magnetic bion" molecule has size r^, 
larger than the 1/L UV-cutoff of the long distance theory (2.11), see [35] for a detailed study. 
Thus, confinement in this theory is generated by this novel type of non-self-dual composite 
topological excitation. Despite being analyzed in the framework of a 3d effective theory valid 
at distances greater than L, the locally 4d nature of the theory is crucial. We note that 
a recent lattice work gave numerical evidence in favor of topologically neutral topological 
molecules [44] . Although this work is done in Yang-Mills theory, such molecules exist in pure 
gauge theory as well. This will be discussed in a forthcoming paper. 

The proliferation of bion-antibion pairs in the vacuum means that the non-perturbative 
zero-temperature partition function of the SU{2) theory with Uf adjoints can be represented 
as a dilute gas of magnetic bions (of charge 2) and anti-bions (of charge -2), which interact 
due to the long range magnetic force: 



Zbi 



N++N- 
bion 



on gas- N.^J 
N±,qa=±2 +• 



N++N- 

n 



L3 



exp 



QaQb 



^ a>b - 



(2.31) 



This is because the "twisted" KK monopoles do not exist in 3d theories: their action, for general expec- 
tation value (Al) — v, is ^ 4^ii) ^ [unit L — >• requires also taking gl 0, with gi — ^ kept fixed. 

Thus the action of KK monopoles is ^ and, in the 3d limit L — > 0, with fixed v and 33, it becomes 

infinite. 
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bion 



bion 




Figure 3: The magnetic bion in the SU{2) theory is a composite of the BPS and KK monopoles, 
which is magnetically charged, but topologically neutral. In Euclidean space, where instantons are 
viewed as "particles", bions should be considered as molecules. The size of a bion is much larger 
than the monopole's, but much (exponentially) smaller than intcr-monopole separation. This allows 
a semi-classical dilute gas approximation for the magnetic bion plasma. 



Here, Ri are the locations of the (anti)bions, the quantity in the exponent is the interaction 
"energy" due to the long-range magnetic field of the bions, and Z^on is a single (anti)bion 
molecular partition function, given by [35]: 



,2 



1 (^, (l+^g2-Loop(^)) 
Zb^on {g{L)) ~ i4_8„ e ^^-Loop' ' , (2.32) 

ffl-Loop W 

where coefficients that play no important role are omitted (the known positive coefficient c is 
also not relevant for our present studies). Using standard methods [29], one then shows that 
the bion partition function leads to a nonzero mass for the dual photon a at zero temperature. 
Thus, the bosonic part of the long-distance effective action (2.11) for QCD(adj) on M^'^ x S 



L 



is modified to: 

2 



S 



2 

1,2 2L\4:TtJ 



{d^a) H ^ cos 2a 

In (2.33), the dual- photon mass is given by [35]: 



a = a + 2^, for SU{2) 
a = a + 4TT, for S0{3) ' 



(2.33) 



The tension 7 of the confining string between test charges is also essentially determined by 



the mass of the dual photon, 7 ~ ^m^}^ 



^It is also useful to give an equivalent representation for the S0(3) theory by using the rescahng a — >■ 2(t 
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Note that if we were to focus just on the bosonic part of (2.6), without considering 
fermions, as in the study of QCD(adj) with massive fermions or deformed Yang-Mihs theory, 
we would have obtained: 



Sh 



^ ^ iS^ \ff^ -^\2 , '^^'^ „™ J / 0- = cr + 27r, for SU{2) 



1,2 2L V47r/ 



{d^a) + — — cos fj 



a = (J + 4^, for 50(3) ' ^'^'^^^ 



In the bosonic SU (2) theory, within the fundamental domain, there is a unique ground state. 
In the 5*0(3) theory, there are two-ground states related by the topological Z2 symmetry. As 
discussed in the Hamiltonian framework in Section 2.2.1, these numbers are doubled when 
one considers the theory with massless fermions. 

Going back to QCD(adj), in the SU{2) theory, recalling (once again) that a = a + 2-k 
(2.19), we find that within the fundamental domain, there are two ground states of (2.33) 
associated with the breaking of the discrete chiral symmetry (2.29) by the expectation value 
of the dual photon: 

1 da k 

Cmin = {0,7r} associated with — / —dz = -, k = 0,1, (2.37) 

2tt dz 2 

where the last equation shows the "winding" number of a, normalized as appropriate for a 
variable of period 27r, between the trivial and A;-th minimum of the potential. The tension 
of the domain wall between the corresponding vacua is proportional to this winding number. 
The domain wall interpolating between the vacua with (a) = and (cr) = vr is associated 
with the spontaneously broken discrete chiral symmetry. 

In the 5*0(3) theory, we recall that a = a + in (2.16). We now see from (2.33) that, 
in addition to the non-anomalous discrete Z2 chiral symmetry (2.29), a— t-cj + vr, the effective 
bosonic V{a) ~ — cos(2(t) potential has an (accidental) Z4 symmetry, which combines the 
discrete chiral symmetry (2.17) with the topological Z2 symmetry cr— )-(T + 27r (2.24). The four 
minima of the bosonic potential in the fundamental domain are now located at: 

I f°° da , k , , , 

fmin = {0, vr, 27r, 37r} associated with — / —dz = -, k = 0,1,2,3, (2.38) 

47r dz 4 

where the "winding" number is now normalized for a field of periodicity 47r. The tension of 
the domain wall interpolating between k = and k = 1 corresponds to the tension of domain 
wall of the spontaneously broken discrete chiral symmetry. The domain wall tension for the 
wall interpolating between k = and A; = 2 is also equal to the string tension Ti for a charge 
probe in the fundamental representation of SU{2). This is because the winding of a between 
the two sides of the domain wall with A; = 2 is the one corresponding to the insertion of a 
Wilson loop of an electric charge 1/2 used to calculate string tensions [29] . 



which makes periodicity into 2n. Then, 5*0(3) action becomes: 



S - 



1 

Trfl,2 2L \2'K 



2 

(O^cr) H — cos 4(7 



a = a + 2Tv, (2.35) 



with four minima within unit-cell. This second representation will also be useful. 
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We end this brief review by noting that the study [8, 9] of this and related Nc > 2 
theories gives the first example of a locally 4d, continuum, and nonsupersymmetric gauge 
theory where confinement can be understood within a controllable analytic framework. In 
particular, the role of the fermions and the importance of non-self-dual topological molecules 
carrying magnetic charge is novel in this regard. 

It is then natural to ask what we can learn about the finite-temperature behavior of the 
theory. 

3. The SU(2)(adj) theory at finite temperature and the Zp XY model 

Before addressing the finite temperature dynamics of the small-L theory, let us recall the 
relevant scales in the zero-temperature problem. The 14^-boson and heavy-fermion Compton 
wavelengths are ~ L, and coincide with the characteristic size of monopole-instantons, r^ . 
The Compton wavelength of the adjoint "Higgs" field — the uneaten "radial" component of 
A4 — is of order ^. The size of the magnetic bions is of order ^ and the typical distance 

OTT 

between bions is Ai?bjo„ ~ Le^a^ , see [35]. Finally, the Compton wavelength (2.34) of the 
dual photon, ^ ~ Le~9^ , is the largest length scale in the problem: 

{m^ ^ rm.) <C <C mon < ARbion < m;;^ . (3.1) 

Needless to say, it is this hierarchy of increasing length scales at small L, illustrated in Fig. 3, 
that makes the problem analytically tractable. 

3.1 The thermal partition function as an electric-magnetic Coulomb gas 

Turning on temperature means that we are considering the Euclidean theory onM^ xgi xS^, 
instead on M'^ ^ ^i- The fermions are now antiperiodic around the thermal circle, but, as 
before, periodic in S^. We will study the behavior of the system with changing P = j, at 
fixed L. At temperatures low compared to the confinement scale (the mass gap), T <^ rrio-, we 
expect the thermodynamics of the theory to be rather trivial, described, for Uf > 1, by the 
thermodynamics of free massless fermions (the components of the fermions along the Cartan 
subalgebra) and a massive a field — the low-T free-field theory limit of the bosonic theory 
(2.33). The thermodynamic potential is S7 = —TlnZ ~ AT + Am^T e~~ , where A is the 
2d volume of the system. As the temperature increases past m„, at first one can still use the 
a free-field theory description, leading to Q given by the high-temperature limit Q ~ AT^ (in 
the limit of neglect of interactions in the dimensionally reduced (2.33)). 

The most interesting behavior occurs as the temperature increases further above nia, in 
the range rrio- ^ T ^ ~ mw-^^ Now recall that at energies below the non-perturbative 
dynamics of the zero-T theory is that of a decoupled (neutral) free fermion and a scalar 
with an exponentially small mass. The scalar sector can be described as a non-relativistic 3d 

^^At even higher temperatures, above mw ~ in the small-L regime we expect to find essentially free-field 
behavior, with the thermodynamic potential scaling as A^^ ■ 
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Figure 4: The scales in the finite-temperature problem. For the interesting temperature range (3.3), 
the bion size is much smaller than inverse temperature which in turn is much smaller than inter-bion 
separation, i.e, < /3 <C Ai?bion- 

Coulomb gas (2.31) of charges ±2 (the magnetic bions). This picture essentially remains in 
the temperature range ma <C T <C ~ mw, albeit with a few important changes, as we now 
discuss. 

First, all fermions — even the massless ones, responsible for the binding of monopoles 
and KK-anti-monopoles into bions — now have a thermal mass of order T. Their Euclidean 
propagators, and thus the attractive potential between the bion constituents, are affected at 
distances larger than t^. Recalling that the bion size is r* = and requiring that the 
behavior of the fermion propagator is only different at distances farger than the bion size 
implies ^ = ^r^- Thus, the fermion-induced attraction between the constituents of 
the bions will be unaffected so long as: 

r< — = ^-^. (3.2) 

The critical behavior that we shall find occurs well within the range (3.2). In fact, a more 
detailed calculation of the attractive potential between bion constituents is possible and one 

2 

finds that near Tc, which, as we show later, is equal to the finite temperature contribution 
to the bion potential near r* is suppressed, relative to the zero temperature potential, by a 
factor of order (r^Tc)^ = {8nf)~^, where r*Tc = l/(8n/) <C 1. Thus, for the dynamics near 
criticality, well within the range (3.2), we can treat bions as pointlike. The ultimate reason 
for this is that the bions are much smaller than the size of the compact "thermal" direction, 
as illustrated on Figure 4. 

Second, recall that the inter-bion separation ARhion at T = is ARhion ^ Le^n ^ 
Thus, for T within the range: 

8'^^ 1 ttL 

Le^ ~ ARuon > 7^ > ~ ^— , (3.3) 

the separation between bions is ^ l/T", (and l/T, from (3.2), is in turn larger than the bion 
size). This means that in the regime (3.3) the bion dynamics, whose T = partition function 



- 23 - 



is (2.31), can be described by the dimensional reduction of (2.31). Thus, the 3d Coulomb 
potential between bions can be replaced by the 2d logarithmic one, where the coefficient 
follows simply from Gauss' law: 

1 ^_2riog|i?|r, (3.4) 

and R denotes now a vector in (the argument T in the logarithm was inserted for dimen- 
sional reasons; it will play no role due to charge neutrality of the 2d gas). Furthermore, the 

9 R ■ (fin ' 

integral over positions of the bions should be now over points in M , replacing — ^-^^ — )• yl^- 
Thus, we obtain for the thermal partition function of the 2d bion gas: 



N++N- N++N- 



Z2,b^on,as= ^ fj^fl U R, e^p 

N±,ga=±l +■ ' j 

where we introduced the bion fugacity: 



SIttLT , - 

2_^qaqb log|i?a - Rb\T 

a>b 



9' 



(3.5) 



2'bion e-2'5o(i+''5) 

■^bion = -JJY ~ l3rpgM-8nf ' (^-6) 

as follows from our comments above (3.5) and from eqn. (2.32); we also note that the Coulomb 
gas partition function in 2d has to obey charge neutrality, A^4. = A^_. We took the sum over 
charges in (3.5) to be over qa = ±1, i.e., the factor of 2 in the bion charge was absorbed in 
the interaction strength. 

Now, recall that the partition function of our theory is Z = ZQZ2d bion gas, where Zq 
describes the free photon fluctuations, now reduced to 2d and Z2d bion gas is from (3.5). Recall 
that the T = partition function Z = ZoZbion gas, when expressed in terms of the a field, 
equals Z = f 'Dae~^ , with S given in (2.33). The dimensional reduction of the T = partition 
function (appropriate in the range of temperatures (3.3)) is equal to ZQZ2dbion gas and is 
precisely the partition function of a XY model in a Villain approximation (see, e.g., [1]).^"^ 
The "vortex-vortex" coupling K can be read off from the coefficient of the logarithm in the 
exponent in (3.5), which describes the interactions of vortices:^* 

K='^. (3,7) 

It is well known from the physics of the XY model that at small values of the spin-spin (and 
vortex-vortex) coupling K < Kc = 4 the spontaneous formation of vortices is entropically 

2 

favored. Note that Kc = 4 is equivalent, from (3.7), to Tc = ^r^. As we will later see, 
this value of Tc is the exact (up to exponentially small corrections) value of the critical 



^^The XY model is defined by a partition function as in eqn. (1.1), but witii y — Q and k — >■ Jf, i.e., 
—PH — X]a;-^=i 2 ^ cos{dx+ii — Sx)- The contribution of the vortices of the 9 field to the partition function 
take exactly the form Z^duon gaa with K defined by (3.7); see Section 3.3 for a derivation in a continuum 
physicist's language. 

^*Note that as opposed to usual discussions of the XY-model, for us small- if corresponds to low-T. 
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temperature for the S[/(2)(adj) theory. At smah K, vortices disorder the system and lead 
to a finite correlation length (mass gap). From the map (3.7), we see that in the S'f7(2)(adj) 

2 

K < A corresponds to the low-temperature confining region T < Tc = At large values 
of the spin-spin coupling K > Kc, the appearance of vortices is suppressed, the dynamics is 
dominated by spin waves (the Zq partition function), and the XY model exhibits algebraic 
long-range order with continuously varying critical exponents and no mass gap. From (3.7), 
this behavior would be attributed to the high-temperature T > Tc regime of the S?7(2)(adj) 
theory and thus be expected to describe the deconfined phase. If (3.5) was all that was 
relevant in this temperature range, it would lead to a high-temperature phase with an infinite 
correlation length — while we expect to have finite correlation length above the deconfinement 
transition due to Debye screening of electric charges. It has already been noted [11], in the 
context of the 3d Polyakov model, that the BKT behavior described above is not the one 
expected of a confinement-deconfinement transition. 



Thus, there is essential physics missing from the description only in terms of the monopole- 
instanton gas (for the Polyakov model), or bion gas (for SU (2)(adj)). The remedy was already 
suggested in [11]: that effects of the heavy W-bosons have to be included in the description 
of the deconfinement transition. While at T ^ niw (see (3.3)) the effects of W bosons are 
Boltzmann suppressed, the effects of monopole-instantons (or bions) are also exponentially 

OTT 

small, with fugacities e ^ . The Boltzmann suppression of M^-bosons {myy = ^ in the 

center-symmetric vacuum) is, at T = Tc, e~ = e ^~a^ = e . Thus, near VF-bosons 
and bions are equally likely to appear in the plasma and we expect that the deconfinement 
transition in S'C/(2)(adj) is driven by competition between the interactions of electrically (VF- 
bosons) and magnetically (bions) charged particles. This is the analytic realization of the 
scenario proposed in [17]. 



Including the effects of electric charges on the bion plasma partition function is relatively 
straightforward. First, note that the neutral particles, such as the radial mode of the "Higgs" 
boson A/^ is not expected to play a role in the deconfining dynamics, despite being lighter than 
the charged VF-bosons. Recall that ^4 is the Wilson line wrapping the spatial, non-thermal 
circle. The heavy charged fermions (of mass equal to / = l,...,nj) contribute 

similar to the VF-bosons, because at ttIo- <C T <C m-vy, the VK-boson, , and bion gas is 
(exponentially) dilute: the thermal de Broglie wavelength (mvyr)~2 is much smaller than 
the typical distance between particles ~ ^Ruon^ and the Fermi statistics is irrelevant. Thus, 
we will further refer to the gas of electrically charged particles as the 'W-boson gas" and will 
account for the A^ contribution via the multiplicity — see also Footnote 20 below. At the T <^ 
mw temperatures of interest, the VF-boson partition function is that of a 2d gas of electrically 
charged non-relativistic particles, which interact via Coulomb forces with themselves (as well 
as with the magnetically charged objects, the bions, as described below). Thus, the W- 
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boson/A^ fugacity is:^^ 

iw = {2nf + 2)'^e-^ . (3.8) 

Another subtlety that we have to discuss is that at the lowest mass level on x 'S>\ 
in the center-symmetric vacuum there are actually two sets VK^-bosons (and fermions). 
This is easiest to see in the D-brane picture [40], which, despite being highly supersymmetric, 
greatly helps in studies of the nonsupersymmetric tree-level spectrum. It can also be seen via 
the usual Kaluza-Klein decomposition: the masses of come from their interaction with 
the Higgs field ^4 (recall Footnote 5), which always enters the Lagrangian as D4, = (?4 + lA^. 
Replacing with ^^^^ (n G Z) and ^4 by its vev diag( J, —j), see eqn. (2.7), we find that 
the masses are proportional to l^!E^±Z!l. Thus states of mass J appear at the n = and 
n = — 1 Kaluza-Klein levels, explaining the presence of two sets of lowest mass bosons. 
These will be the only Vl^-boson states whose contribution we will keep. Similar to accounting 
for the A^ contribution, at T ^ mw, the contribution of these states to the grand partition 
function can be accounted for by doubling the fugacity (3.8).^*^ 

Because the Vl^-bosons carry electric charges, two Vl^-bosons with electric charges qa and 
Qh (equal to ±1) located at Ra and Rh interact via the logarithmic potential:^^ 

Vwab = - ^ J rj. (laQb ^'a.\Ra - Rb\T . (3.9) 

In addition to the 2d Coulomb interaction (3.9) between VF-bosons, there is an Aharonov- 
Bohm phase due to the presence of magnetic charges (the bions) in the system. The interaction 
between the a-th bion located at the origin in and the 6-th static VF-boson of electric charge 
Qh, located at (in M^) is given by the integral over the VF-boson worldline (i.e., along xq): 

l/T 

dxoAO/„„(xo,X6) , (3.10) 

where ^^■^^(a^O)^) is the time component of the gauge field of the bion. Note that the 
exponential of i times (3.10) contributes a phase factor in the path integral also in Euclidean 
space and that in our normalization of the fields no factors of g appear. Next, we note that 



^^The prefactor follows from integrating over the momenta in the non relativistic partition function, equal to 

the product of factors j (27!)'^^ e~ t~ 2nir +^(^) for all particles. The 2n/ factor accounts for the multiplicity 
of charged fermions. The various terms contributing to the interaction V{x) are given in eqns. (3.9) and (3.11). 

^"Perhaps a comment on this is necessary. We are treating the two kinds of bosons (and the A* fermions) , 
which have the same charges (and, to the order of our calculation, the same fugacities), as indistinguishable. 
One can show, via the sine-Gordon representation of a 2d Coulomb gas partition function, using the fact 
that only overall charge-neutral configurations contribute to the 2d partition function, that, indeed, the grand 
partition function of a gas of two kinds of same-charge particles with fugacities ^1 and ^2 is equal to that of a 
gas of one kind of charged particles with fugacity + £,2- 

quick way to obtain this formula is to recall that unit electric charges appear as unit winding number 
vortices of the dual photon field in (2.1); then (3.9) is just the interaction energy of two vortices of unit 
winding. Equally quickly, since the interaction energy of two static W bosons in R^'^ is 
W^s propagate in a (Euclidean) time interval the corresponding action is precisely (3.9). 
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on ]R2 X Sj^ the field that the ty-boson experiences is that of an infinite chain of bions (ah 
located at the origin in M?) and spaced equidistantly along the xq axis. Thus, we can replace 
(3.10) with: 



Qb 



/ dxQ ^bion{^0 + 7^.Xb) = qb dxo^^iKxcXb) =4gfee(ffe) , (3.11) 

Jo ^ J-oo 



where Q{xa) can be taken to be the angle between the vector Xa connecting the monopole 
and the VK-boson and the positive X2 axis. To obtain the last equality above, we arranged 
the Dirac strings of all monopoles to be along the negative X2 axis and used the gauge field 
of a monopole of charge Qa (recall Qa = 2 for bions). In usual polar coordinates it is given 
by ^ = Qa(l — cos9)d(j), when the Dirac string is along the negative z-axis. This expression 
for A is then adapted to our coordinates by replacing x — )• xq, ?/ — )• xi, z — )• X2 (thus the 
strings are now along the negative-X2 axis). Using it to calculate the integral yields (3.11); 
charge neutrality was also used to drop some irrelevant constant contributions. The heavy 
(mass mw) fermions X± have the same Aharonov-Bohm interaction with the bions as the W 
bosons and their spin-orbit interaction with the bion magnetic field is mass-suppressed. 

Thus, including the VF-bosons, in addition to the bions, we arrive at the following par- 
tition function (note that, as in (3.5), we absorb the factor of 2 from the bion charge in the 
prefactor of the bion-bion interaction) : 

^bion+W 



Nb±,qa=±l Nw±,qA=±l 



X exp 



2 



^'^^^'^ q^q^ In \Ra - Rh\ + i^^jjf IaQb In \Ra - + 4f ^ qsqaQiRs - Ra) 

^ a>b ^ A>B a,B 



The meaning of the various terms in the partition function above have already been explained. 
Eqn. (3.12) is the nontrivial (i.e., interacting) part of the partition function of the 5[/(2)(adj) 



theory on M^'^ ^ gi , for T in the range j^e < T < as in(3.3). 

An important property of Z^ion+w is that it is invariant under electric-magnetic duality. 
This follows immediately by noticing that exchanging: 



S.bion '■' '' 2^^^/, 

327rLr c/2 



2ttLT 



(3.13) 



2 

in (3.12) gives rise to an equivalent partition function. The self-dual point T* = occurs 
exactly at the critical temperature of the bion-only gas. This temperature also happens to 
be the BKT temperature of the VF-only gas (where the only interactions would be given by 
(3.9)). This strongly suggests that the deconfinement transition in S'C/(2)(adj) indeed occurs 
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The partition function Zfjio^+w^ when defined on the lattice, also has a Krammers- 
Wannier type duality, analogous to (3.13) (exchanging electric and magnetic charges living on 
dual lattices) , and is known [6] to be equivalent to the lattice XY- model with a Z4 preserving 
perturbation, defined in (1.1). We will refer the reader to the quoted literature for the lattice 

duality; instead, will present a (somewhat shorter but helpful) continuum version later in 

2 

Section 3.3. There, we will also establish that k = 2^lt ^ claimed in Section 1.1. 

Before we continue, we note that if one studies the thermal physics not of QCD(adj) but 
of deformed Yang-Mills theory on M^'^ x [12], one finds, instead, a lattice XY- model with 
a Z2 preserving perturbation. If the 3d minimally supersymmetric Georgi-Glashow model 
is studied at finite temperature, a partition function similar to (3.12) is obtained [46]. The 
nature of the composite topological excitations in the A/" = 1 3d theory, analogous to magnetic 
bions, was only elucidated recently in [47]. 

3.2 The renormalization group equations for the magnetic bion/W-boson plasma 
and the approach to Tc 

The behavior of the magnetic-bion/W-boson Coulomb gas (3.12) can be studied by various 
means. One way would be to do Monte-Carlo simulations, using its representation as a 
lattice XY-model with a Z4-preserving perturbation, discussed in Section 3.3. Another way 
is to study the perturbative renormalization group equations and look for a fixed point where 
the correlation length diverges. It is interpreted as the critical point corresponding to a 
continuous phase transition. Showing that this point exists for the Coulomb gas (3.12) is the 
subject of this section. 

To lowest order in the fugacities £,uon and 2^w, these equations are derived in a physically 
intuitive way in Appendix A. There exists much literature on the subject, see, e.g. [1,3,48-51], 
but we find that to the order we are working the derivation given in the Appendix is the 
simplest we are aware of. To describe the renormalization group equations (RGEs), we 
introduce the dimensionless variables: 



where a is the "lattice spacing" (inverse UV cutoff, ~ L) used to define dimensionless fugac- 
ities. In terms of these variables, the duality relations (3.13) become: 



To leading order in the fugacities, the RGEs (A. 4, A. 5, A. 14), given here for the Z4 model are: 



y = 2iw a^, y = ibion a^, k 




(3.14) 



16 

y ^y, K o — 

K 



(3.15) 




(3.16) 
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where y = and the course graining parameter is a — )• e a with 6 > 0; thus increasing b 
corresponds to flow to the IR. The first equation above describes the screening of the Vl^-boson 
Coulomb interaction of strength k by electric charges (the —K^y^ term on the r.h.s.) and its 
anti-screening by magnetic bions (the term) . The origin of the antiscreening term is in the 
Aharonov-Bohm interaction between VF-bosons and bions, see discussion above eqn. (A. 14) 
in the Appendix. The last two equations in (3.16) describe the change of the fugacities of 
magnetic bions (A. 4) and VF-bosons (A. 5) upon coarse-graining. 

The RGEs (3.16) are invariant under the electric-magnetic duality (3.15). Thus, the two 
terms responsible for the "running" of the magnetic coupling — : 



anti-screening by VF-bosons. Since the duality is an exact symmetry of the partition function 
(3.12), it is also a symmetry of the exact RGEs. 

At the self-dual point of (3.15) we have k = A, y = y. Thus, for the Zp=4 case of interest 
(3.12), this point corresponds to a line of fixed points of the leading-order renormalization 
group equations (3.16). The duality (3.15) maps high to low temperatures and, as usual, 
suggests that if there is a unique transition, it should occur at the self dual point. Thus, in the 

2 

Z4 theory, the deconfinement transition is expected to occur at the self dual point Tc ~ g|r^ , up 
to small corrections (the precise equation determining is given in (B.8)). The exponential 
accuracy of this determination of Tc as well as the flow to the critical point are described in 
detail in Appendix B. The important point to make here is that the perturbative expansion 
in small fugacities in the S'C/(2)(adj) theory is sufficient to reliably study the approach to the 
critical point, because the fixed line extends to zero fugacities. This is in contrast with the 
Z2 model studied in [11, 12] where the small- fugacities approximation breaks down, but is 
similar to [46] , where it is mentioned but an analysis of the RG flow is not discussed in detail. 

The study of the RGEs near the fixed line is considered in detail in Appendix B, where 
we show that the correlation length diverges as:^^ 



where yo and yo are the VF-boson and bion fugacities taken at the UV-cutoff scale. Hence, 



the free energy F ~ ^"^ ~ [T — Td ^^^^0^0 ig continuous at the transition along with its 
derivatives to a large (but finite, unlike the BKT transition) order. 

The critical exponent 77, measuring the decay of the correlations at T = Tc, also depends 
on the values of the fugacities on the line of fixed points. It has been shown [4], using the 
dual-sine-Gordon description of the p = 4 Coulomb gas, that the theory at the self-dual point 
is equivalent to that of a free scalar field (2d c = 1 conformal field theory), with critical 
exponents that depend on the fixed-line value oi y = y (an earlier perturbative calculation 

As far as we can tell, this result is new; however, it may exist somewhere in the condensed-matter literature. 




(3.17) 



written in the dual frame, can be understood as being due to its screening by bions and 




(3.18) 
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Figure 5: A two-dimensional (in the y—n plane) cross-section of the RG flow of the Z4 clock model for 
small fugacities, near the critical coupling k^, = 4. A two-dimensional flow is appropriate, since yy is 
approximately RG invariant near Kc, see (B.3). When the temperature is tuned to the theory flows 
along one of the two critical trajectories (shown in thick lines) to a position on the fixed line of (3.f6): 
y = y = y^/An, k = 4. determined by the UV values of the fugacities (c' = IGTT^yoyo; in the plot, 
we have chosen \/d = 0.2). The equation determining and the critical trajectory corresponding 
to 47ryo > (the one in the upper right quadrangle) are found in Appendix B. Note that when 
one fugacity becomes small the other one grows and that both fugacities remain small only on the 
critical trajectory. The amount of RG "time" that trajectories stay close to the critical one diverges 
asT Tc and is used to find the divergence of the correlation length given in (3.18), similar to the 
analysis of the BKT transition. 

of the conformal anomaly, also indicating c = 1, is given in [49]). We will not need to make 
use of this result, as the fugacities in the small-L theory are already small; hence aX T = 
the decay of correlations is governed by the BKT exponents with = 4, with negligible 
corrections due to the nonzero y- 

3.3 Dual descriptions and symmetry realizations above and below Tc 

In this section, we study the duality between the Coulomb gas (3.12) and XY models with 
symmetry-breaking perturbations. 

3.3.1 Zp XY models as Coulomb gases 

The lattice XY-model with a Zp preserving perturbation is defined by a partition function 
similar to (1.1), but with an "external field," breaking the continuous U{1) symmetry 9^ — )• 
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9x + c to a Zp subgroup, 9^ 9x + y: 



2n 

ZxYZ, = I lld9i e ' , (3.19) 

* 

where g is the strength of the C/(l) — )• Zp perturbation (the lattice spacing a is set to unity). 
We have denoted the couphngs by K and g (instead of k and y as in (1.1)), since we are going 
to use the duahty of a Coulomb gas to (3.19) in several different ways. 

We will not derive the duality of (3.19) to (3.12) in the most rigorous way, as a detailed 
derivation on the lattice can be found in the literature [3,6]. Instead, we will give a more 
qualitative continuum presentation, see [22,52]. While less rigorous than [3,6], it is perhaps 
more intuitive to the continuum physicist. 

To begin, note that the naive continuum limit of (3.19), obtained by expanding the cosine 
and replacing it with the 9 kinetic term (shown below) , is the Euclidean theory of a compact 
scalar 9 = 9 + 2tt with Euclidean action and partition function: 

S[K;G,p]= [ ^{di9f -2G cos p9 , Z[K-G,p;H,w]=[v9e-^^^'^^'P^, (3.20) 
Jr2 47r J 

where 2G = ^ and a is the lattice spacing. The meaning of most the arguments of Z is clear, 
except for w and H. These will be explained in detail below (but let us mention that w is 
the (minimum) winding number of vortices that appear with nonzero fugacity H). To argue 
that Z[K; G,p; H,w] maps to an e-m Coulomb gas, begin by expanding the interaction term 
as follows: 

fe>o ■ V-^ ^ / 

Then, noting that performing a Gaussian integration over 9 will make only equal numbers of 
e*^ and e~*^ contribute (this imposes charge neutrality on the resulting Coulomb gas), rewrite 
(3.21) as: 

^fd^.Gco.pe^Y^ E / rf'^Ae^^'^'^"^^ ■ (3.22) 

LI '"-•) J 



Ea '3A=0 



(n!)2 11 

n>OqA=±l A=l ■ 

We then insert (3.22) in (3.20), interchange the orders of the sum over n and the path integral 
and perform the Gaussian integral over 9. Clearly, every term on the r.h.s. of (3.22) is now 
a source of 2n "electric" charges pqA at positions xa- However, the most general solution for 
9{x) also includes a set of an arbitrary number of vortices, allowed due to the periodicity of 
9. Take these to have winding numbers qa and be located at positions Xa-'^^ Note that the 
gas of vortices must also obey a neutrality condition Y2a 9a ~ ^' otherwise the action will 
diverge. Thus, the general solution for 9{x) in a sector with a given number of vortices and 
charges is: 

9{x)ciass = --^^QA^T^ \x - xa\ + ^qa&{x - Xa) + 9o{x) , (3.23) 
A a 

^^Including vortices requires a UV definition, provided, e.g., by the lattice model. 
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where 9o{x) represents periodic spin- wave fluctuations and Q{x) is the polar angle between 
X and, say, the positive x axis. We note that as far as vortices are concerned, we are free 
to consider vortices of arbitrary winding number, with fugacities that can depend on the 
winding. In what follows we will only sum over vortices with \qa\ = w. Note also that while 
the vortex fugacities do not explicitly appear in (3.20), they are part of the definition of the 
theory, and are required by the 2d compact scalar duality [52]. 

The final step is to substitute the classical solution in (3.20) and calculate the action 
(omitting all self energies which are, of course, absorbed in the fugacity after renormalization). 
The spin- wave fluctuations decouple from the vortices and the charges, while the latter only 
interact via an Aharonov-Bohm type interaction, which arises due to the phase factors in the 
expansion of the cos p6 term. The flnal result, including also a sum over the arbitrary number 
of vortices of fugacity H, gives the partition function of the compact scalar theory (3.20) in 
the form: 



Z[K;G,p;H,w] = Zox (3.24) 

/^2n Tj2m f 



2 

r , ^ ^ n ^A^j^ a>b A,b 



m,ga=±'w n,qA=±i 

where Zq is the partition function of the massless scalar Oq representing the spin waves. The 
electric and magnetic (winding number) neutrality is understood in (3.24) and the sums over 
n, m are from to oo. The integral over the position of the charges and vortices is denoted 
by J [(i^("'"'""^V] . Finally, as written, the partition function only involves a sum over one set of 
magnetic vortices, but the generalization to many is trivial. 

The summary of this section is that we found that the compact-^ theory (3.20), related to 
the XY-model (3.19), has a representation in terms of an e-m Coulomb gas (3.24). In what 
follows, we shall call (whenever convenient) the charge-p particles labeled by qa, xa "electric" 
and the winding- number w vortices labeled by qa^fa — "magnetic". 

3.3.2 Zp XY-models as duals of the W-boson/bion gas 

Now we can map the bion-/VF-boson gas partition function of the SU{2) QCD(adj) theory 
to the Coulomb gas (3.24) of the Z4 XY model. To facilitate the comparison, we reproduce 
here the interaction energy ("— /3//") of the VF/bion gas from (3.12): 

— E Q'agfelnrafe + K ^ qAqE'^'^^rAB+'^i^^qsqaQirBa) , ga,A = ±1. (3.25) 

a>b{bions) A>B{W bosons) a,B 

There are two ways to map (3.24) to the W^/bion Coulomb gas (3.12) that are inequivalent 
under e-m (vortex-charge) duality. The flrst matches to the SO{3)=SU{2)/Z2 gauge theory 
and the other matches to the SU{2) theory. We discuss both in turn: 

1. First, we can identify the charge-1 particles with the VF-bosons of (3.12), and the charge- 
p particles with the magnetic bions. Upon comparing to the bion partition function 
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(3.12), or the "energy" displayed in (3.25), we find the map between (3.20) and the 
QCD(adj) theory: 

p = /^^ yj = l^ K = K= G = Cbiorn H = 2Cw. (3.26) 

where k is as defined in (3.14). Thus the dual description is the theory defined by (3.20) 

with Z[K\CuonA'iKwA\- 



s, = 



[ -^{diOf -2Cbion COS ie, z[K-XbionA]Kw,i]= I ve e~^^ . (3.27) 

Jr2 47r J 

Note that the bion-induced potential given in (2.33) for SO{'i) is Att periodic. On the 
other hand, our angular spin variable is 27r periodic. It is, therefore, more convenient 
to compare to the equivalent alternative form for 5*0(3) given in (2.35), which agrees 
with (3.27). 

Since unit-winding vortices are identified with VF-bosons and charge-4 particles — with 
magnetic bions, which carry magnetic charge Qm = 2, the description by means of Si 
allows us to also introduce particles of integer charges < 4 into the system. In the 
language of the original gauge theory, this means that we could introduce monopoles 
of charges as low as Qm = \ — the minimal charge allowed by Dirac quantization in the 
50(3) theory, equal to one-quarter of the bion magnetic charge. On the other hand, 
since the M^-bosons are the unit- winding vortices and since there are no ^ -winding 
vortices — as these would violate Dirac quantization — we can not introduce dynamicaP^ 
particles oi Qe = \ into the system. Thus, the description (3.26, 3.27) is appropriate 
for the 50(3) theory. 

The e-m dual of (3.26, 3.27) is obtained by interchanging p and w, i.e., by letting p = 1, 
w = A, and taking the VF-bosons to be the "electric" charge-p = 1 particles in (3.24), 
and the bions — the w = 4-winding "magnetic" particles. To get the interactions in 
(3.12) reproduced correctly, this requires taking: 

p = l, w = A, K=^ = G = 2Cw, H = Cuon . (3.28) 

We note that (3.28) is simply the transformation of (3.26) under the e-m duality (3.13) 
and is also the usual T-duality of S"*^ sigma models (compact-^ theory) in 2d [52]. T- 
duality interchanges K -H- l/if, i.e., inverts the "radius" of the compact field 9, and 
also interchanges vortices and charges (or "electric" and "magnetic" particles, in the 
Coulomb gas language). Overall, the electric-magnetic duality amounts to: 

p^w, H^G. (3.29) 

K 



■^"^ However, the theory can still be probed by external charge- 1/2 "electrons" (despite the fact that they 
would violate Dirac quantization if they were dynamical). In the 2d Coulomb gas picture, this corresponds 
to the insertion of, say, two external charge-1/2 "magnetic" particles in the plasma and studying how their 
interaction is affected by the dynamics of the particles in the plasma as the temperature changes. 
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This dual description is thus the one defined by (3.20) with 2(w, 1; Cuon-, 4]: 

Si = [ -^{d^e)^ -4.Cwcose, z[-;2Cw,i;CbionA]= j ve e-'^^ . (3.30) 

J^2 Attk k J 

Both (3.27) and (3.30), the model with action Si and its dual S'l-model, describe the 
same 50(3) theory, but are weakly coupled in difi^erent regimes. 

For later use, note that we use to denote the spin variable in (3.27), such that e*^^ 
creates magnetic charge in the language of the original gauge theory — minimal charge 
monopoles for q = 1, 't Hooft-Polyakov monopoles for q = 2, and magnetic bions for 
q = 4. In contrast, we used 9 when describing the T-dual spin model in (3.30), such 
that e^*^^ creates electric charges in the original gauge theory, which are H^-bosons for 
q = 1 (as this is the 5*0(3) theory, "electrons" are not allowed as dynamical objects). 

2. The second interesting possibility, which is not related by 2d e-m duality to the pre- 
ceding one and relevant to our later discussion, is to take the p = 2 charges in (3.24) 
to correspond to the magnetic bions. To preserve the Aharonov-Bohm interaction, we 
must now take the unit winding {w = 1) charges in (3.24) to have zero fugacity, and 
instead sum over charges with w = 2 (assigning them fugacity H). This restriction 
due to the Aharonov-Bohm interaction is protected by a Z2 symmetry of the theory, 
and consequently, w = 1 vortex configurations are forbidden. The "magnetic" charge-2 
particles of (3.24) will be now identified with the Vl^-bosons. We then have another 
representation of the bion/VK gas (3.12) in terms of: 

2 

p = 2, w = 2, K = j = G = iuon ,H = 2Cw . (3.31) 

This dual to the bion/1^ gas is the theory (3.20) with (bion, 2; 2(y[/, 2]: 

•^2= / r^idiOf -2^Mon COS 26 , Z[^;^bion,2;2^w,2]= [vee-^^ . (3.32) 
Jr2 16vr 4 J 

In the description (3.31, 3.32), the VF-bosons are winding number-2 vortices and the 
magnetic bions are the charge-2 particles. The theory permits the introduction of 
Qe = \ probe "electrons" (which would be the winding number- 1 vortices) as well as of 
Qm = 1 't Hooft-Polyakov monopoles as dynamical objects (the latter can be introduced 
by adding Majorana masses for the gauginos to lift the fermion zero modes). Thus, the 
description in terms of 5*2 is appropriate to the SU{2) theory. Indeed, (3.32) is the 
dimensional reduction of the bion induced potential (2.33) for SU{2) down to M^. 

To construct the e-m dual, we can take the p = 2 charges in (3.24) to correspond to the 
M^-bosons and identify the w = 2 vortices of (3.24) with the magnetic bions. Then we 
have another representation of the bion/VF gas (3.12) in terms of: 

p = 2, w = 2, K = - = G = 2Cw, H = iuon ■ (3.33) 



- 34 - 



Thus, now the theory dual to the bion/VF gas is (3.20) with 2(w, 2; Cbiom 2]: 



S2 




{diOf - Aiw COS 29, Z[-- 2iw. 2; iuon, 2] 



/ 



(3.34) 



This theory allows the introduction of dynamical "electric" charge- 1 particles (the "elec- 
trons" of the original gauge theory with Qe = |) as well as "magnetic" charge- 1 particles 
(which in the gauge theory correspond to 't Hooft-Polyakov monopoles with Qm = 1) in 
the plasma. However, in (3.34), we can not introduce dynamical magnetic monopoles 
corresponding to Qm = \- Thus, the description (3.33, 3.34) is appropriate for the 
SU (2) gauge theory. 

We stress again that the difference between the two descriptions above lies in the allowed 
magnetic and electric charges of dynamical objects that can be added to the theory without 
violating Dirac quantization. We caution the reader against thinking that Si/ Si and S2/S2 
map to different — from each other and from (3.12) — lattice models. This is not so, provided 
the unit- winding fugacity in (3.33) is strictly zero (the models differ in the charges of allowed 
dynamical objects not present in the VF-/bion-gas (3.12)). Their equivalence as Coulomb 
gases follows from our discussion and, e.g., the lattice representation [6].^^ 

In each case — either SO{3) or SU{2) — we used e*''^ to describe the creation of magneti- 
cally charged objects in the gauge theory, and e**^^ — to create electrically-charged gauge-theory 
objects. Various correlators of these exponentials probe the phase structure of the theory and 
will be discussed in the next section. 

To summarize this section, we have established two (up to e-m duality) ways, appropriate 
to SO{2>) or SU{2) gauge theories, to map the 5C/(2)(adj) partition function (3.12) to the 
lattice XY-model with a Zp-preserving perturbation (3.19). In each case, one side of the 
duality is a thermal 4d nonabelian gauge theory with massless fermions, compactified on a 
small circle, while the other side is a well-studied lattice model. The dynamics of the 4d 
QCD(adj) gauge theories is complicated and not well- understood. Numerical Monte-Carlo 
simulations of theories with massless fermions are expensive and time consuming. On the 
other hand, in the simplifying regime of a small-L compactification, theories with adjoint 
Weyl fermions have yielded to controlled theoretical analysis. We have taken this analysis 
one step further by having shown that the thermal dynamics of the theory in this regime can 
be studied via a lattice spin model. 

3.3.3 Topological, center, and chiral symmetry realization 

The dual representations described in the previous section can be used to study the realization 
of the symmetries as a function of the temperature. 

■^^One can show that a Coulomb gas of electric particles of charges quantized in units of of a positive integer 
ke and magnetic particles with charges quantized in units of a positive integer A:,„, with Coulomb interaction 
strength governed by k, is equivalent to a Coulomb gas of electric particles quantized in units of fc^ = 1 and 
magnetic particles quantized in units of fc^ = kekm and an interaction strength by showing that 

they map to equivalent lattice models; that this is so should also be intuitively clear from the fact that they 
have identical values of the "energy" —jSH. 
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The SO{3) theory: The dual descriptions are given by (3.27) and its T-dual (3.30). In the 
description by 5i of (3.27), fluctuations of 9 are small at large values of k, which from (3.26) 

2 

corresponds to low temperature T < Tc = £^ . Then, the state of the theory is determined by 
minimizing the classical potential ~ ^bion cos AO and corresponds to the spontaneous breaking 
of the Z4 symmetry (recall from the T = discussion near eqn. (2.38) that this corresponds 
to spontaneously broken topological Z2 and chiral symmetries). 

On the other hand, the e-m (or T-) dual model described by Si (3.30) becomes weakly 
coupled at large values of i.e. at high temperature T > Tc. The classical potential is cos0, 
with a unique ground state and no broken symmetry in the high-temperature regime of the 
50(3) theory 

The symmetry realization can also be inferred from studying the correlators of the order 



parameter e**^ for Z2 and the 'L'^ chiral symmetry; recall that e**^ is an operator creating a 
minimal charge monopole in the plasma (which is really a 4d minimal charge 't Hooft loop 
winding around the S^). Below, we show the e*''^, q = 1,2,3,4, correlators, as well as those 
of the Qe = 1 Polyakov loop (e*^, representing a VF-boson worldline winding around the S^). 
The Qe = 1 Polyakov loop is not associated with any symmetry, and its correlator is shown 
for completeness.^^ Recall that in the SO{3) theory, there are no center symmetry, Qe = ^ 
charge, and fundamental Polyakov loop. 



SO (3) theory 


T <Tc 


T = Tc 




Qm= 't Hooft loop : (e*^^^^'^e-*3^W) 1^1^^ 


1 


1 




Qe = l Polyakov loop : (e^^»e-*^'(°)) .|^|^^„ 




1 

|2:|4(l+0(y)) 


1 


7j^^^* chiral and topological symmetry 


broken 


c = 1 CFT 


unbroken 


magnetic charges 


m— free 


m— free 


m— confined 


electric charges 


e— confined 


e— free 


e— free 



(3.35) 

The table in (3.35) shows the spontaneous breaking of (and Z2 chiral) at low temper- 
atures, the almost-BKT scaling at Tc, and the confinement of monopoles in the deconfined 
phase. Both the string tension ct(T) and the dual string tension dq{T) are proportional to 
the spin-model mass gap (the inverse correlation length ^~^)., and both vanish as T — t- Tc- 



T 



T 



\T-Tr 



(3.36) 



The value of the exponent follows from our RG analysis, see (3.18) and Appendix B.2; recall 
that the RG analysis is valid on both sides of Tc, owing to the high-T/low-T duality. 

The SU (2) theory: We begin with the 52 theory (3.32), and observe that it is weakly-coupled 
at large k, i.e., at low temperature T < Tc. Thus, semi-classical analysis is appropriate. It 



For the Qe = 1 Polyakov loop correlator the "|a;| — >■ 00" limit should be understood in the sense that while 
|a;| is large, it is still sufficiently small so that the Polyakov loop is not screened by pairs "popping out" 
of the vacuum and leading to string breaking. 
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shows that the potential cos 20 has two minima, indicating that the Z2 chiral symmetry is 
broken in the low-temperature phase. The order parameter is e*^, creating a 't Hooft-Polyakov 
monopole (a 't Hooft loop of a Qm = 1 monopole winding around S}^; it is charged under 
the chiral symmetry due to the intertwining of the topological shift symmetry with chiral 
symmetry, eqn. (2.29)).^'^ On the other hand, the description by means of §2, (3.34), T-dual 
to (3.32), is semi-classical at small k, high temperature (T > Tc), and shows that there are two 
vacua of the cos 29 potential. The corresponding order parameter is the "electron" creation 
operator e*^^^-* — the Polyakov loop of a Qe = 5 "electron" with a worldline winding around the 
thermal circle — whose expectation value indicates that the center symmetry is broken 
at high temperature. We can now summarize the relevant correlators — of 't Hooft-Polyakov 
monopole operators, e*^, and of "electrons", e*^ — in the SU{2) theory as follows: 



SU{2) theory 


T<T, 


T = T, 


T>T, 


Qm = l 't Hooft loop : {ei<^i-)e-^m) ^^^^^ 


1 


1 

\x\l + 0{y) 




Qe = I Polyakov loop : (e^'(")e-*^'(o))||^|^^ 




1 

\x\^+0{y) 


1 


Z2 center symmetry 


unbroken 


c = 1 CFT 


broken 


1^2^ discrete chiral symmetry 


broken 




unbroken 


magnetic charges 


m— free 


m— free 


m— confined 


electric charges 


e— confined 


e— free 


e— free 



(3.37) 

which is the SU{2) equivalent of the SO{3) table (3.35). 

4. On the deconfinement transition in Nc > 2 QCD(adj) 

In this section, we study the deconfinement phase transition for SU (Nc) QCD like theories 
with Uf adjoint fermions compactified on M^'^ x S]^. The SU{Nc) gauge symmetry is broken 
by the Higgs mechanism down to C/(l)^'=~^, as described in Section 2.1. 

4.1 The Coulomb gas for Nc > 2 and Nc = 3 electric- magnetic duality 

The non-perturbative dynamics is also similar to the one described in Section 2.2 for SU{2), 
but one has to take into account the fact that there are now different kinds of non-perturbative 
monopole-instanton configurations. The monopole-instantons are now labeled by the affine 
roots of the SU{Nc) Lie algebra. Let the simple roots be a,, i = 1, . . . A'^c — 1. A given root 
di can be described as an A'^c — 1 dimensional vector. The affine root a^^ is defined as: 

= -ai - 02 - as - • • • - oaTc-i • (4-1) 

^^We warn against identifying e'*-correlators in the SU (2) and SO{3) dual models: a Qm = 1 't Hooft- 
Polyakov monopole is created by e** in the SU (2) dual and by e'^" in the SO(3) dual; similar care should be 
exercised in comparing the e'^ correlators. 
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The extended (or affine) root system A^^ = {a^, a^, . . . cS^''}, using normalization Trtatb = 
^ for the generator of the Lie algebra as stated in Section 2.1, obeys: 

Si ■ OLj = 6ij - ^ij+i - ^Sij-i, with i,j = I... Nc{Nc + 1 = 1,0 = Nc) . (4.2) 

The self-dual BPS and KK monopoles in SU{Nc) are labeled by the extended roots — or 
equivalently, their magnetic charges under the unbroken f/(l)^'^~^. The magnetic charge of 
each monopole is found by integrating its U{1)^''~^ magnetic field over spatial infinity in M'^: 

f An ■ ^ 

/ dS • F = — a* = AttQm- , for the type-i (= 1 . . . Nc) monopole . (4.3) 
Js^ 9 



,2 



In the center-symmetric vacuum, the action of each monopole (4.3) is e "^^ = g g'^Nc . 
Recall that the charges of the VF-bosons are Qwi = thus, 

( 1 i = j 

QMrQw, = \-l ^ = i±i • (4-4) 

otherwise 

The Dirac quantization condition is saturated for the nearest neighbor charges on the Dynkin 
diagram, and is twice (as in SU{2)) the minimal bound for charges on the same site on Dynkin 
diagram. 

Magnetic bions, the topological excitations responsible for the confinement, also exist for 
SU{Nc) with massless adjoint fermions. Bions are composed of a monopole with charges ai 
and anti- monopole with charges — aj_i. The monopole constituents of the bions repel because 
of Coulomb interaction and attract because of fermion zero- mode exchange; this ensures the 
stability of the bions, as in the SU{2) case. Clearly, there are also bions, labeled by the 
position of (say) the monopole constituent on the extended Dynkin diagram. Every bion 
carries a charge under the unbroken ^7(1)^'="^ gauge group, given by: 

Qi = ai- OLi-x , i = l...Nc . (4.5) 

The (magnetic) Coulomb potential energy between two bions (or antibions) carrying charges 
Qi, and Qj, and located at and Rj, is given by: 

9 I Ri — Rj I 

The products of the bion charge vectors can be evaluated, using (4.5) and (4.2), yielding: 

Qi ■ Qj = 3(5ij - 25ij+i - 2(5jj_i + -Sij+2 + 2^i,j"^ lor Nc > 5 , 

Qi ■ Qj = 35ij - 25ij+i - 25ij-i + 6ij+2 for A^^ = 4 , (4.7) 

^ ^ 3 3 

Qi ■ Qj = 3(5ij - -Sij+i - IotNc = 3 , 

Qi ■ Qj = 46ij - 4Ji j+i for = 2 . 
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Note the fact (of later importance) that, for SU{3), we have: 



Qi -Qj = 3 Si - dj . 



(4.8) 



At temperatures mo- <C T <C mw, where mo- is the dual photon mass, the inter- 
bion potential should be replaced with a logarithmic one, as was done for the SU{2) case, see 
(3.4). Hence, the pairwise interaction between a collection of bions (with charge Qi) and 
fii antibions (with charge —Qi) is given by:^^ 



SttLT 



bions 



9^ 



^ ^ Qi ' Qj 



ni,nj 



a>b 
rii fij 



a>b 



(4.9) 



Next, we turn to the W bosons. In Dynkin space we have N(. W^-bosons at the lowest Kaluza- 
Klein level, for reasons similar to the ones explained for the SU{2) case. Each P^-boson 
carries a charge given by gcti , and each has mass M\y = . The would-be 3d interaction 
between two M/^-bosons is: 



Vw 



g aj ■ aj 
'^T^\Ri-Ri 



(4.10) 



At temperature T, the interaction between m^ (with charge di) and mj W (with charge 
—di) bosons is: 



Sw 



2'kLT 



^di 

i>j 


dj 


mi,mj 

E 

A>B 






mi,rhj 

-E 

A,B 



mi,mj 



+ ln(T\R 



•-A ^B' 



A>B 



(4.11) 



The Aharonov-Bohm phase interaction between a Vl^-boson (with charge dj at yj), and a 
bion (with charge Qi at Xi) reads, similar to (3.11): 



2Qi ■ dj@{xi - ijj), 



(4.12) 



where &{xi — yj) is the angle between the vector connecting the monopole and the VK-boson 
and a chosen spacial direction. Hence, the interaction between nij and fhj VF-bosons (with 
charge di and —di respectively), and rii and fii bions (with charge Qi and —Qi respectively) 



Note that for Nc = 2, there is only one bion, hence there is no sum over Nc diflerent kinds of bions. 
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is given by: 

'S'vFjbions 



a,B 



ni,mj 



rii ,mj 



a,B 



B,a 



(4.13) 



Collecting everything, we arrive to the 2d Coulomb gas representation of the partition function 
of SU{Nc) QCD(adj): 



Zhion+W = ^ 



,w 



0+ 0— • 



"w+"w~ 



n n 



X exp [— 5bion — Sw — Sw,Uons] , 



(4.14) 



where ^bion and S,w are the bion and W fugacities ^bion ~ / ^4~an,f , ~ rnyjTe 
For further use, let us rewrite (4.14) as follows: 



{Nl^>0, i>0, qa=±l) {N^±>0, i>0, gA=±l) 



(^) 



11 m+ »n— e+ rn 



X exp 



qAQBai ■ dj In — ^ H ^1 ga^bQi • Qj In — 

i>j A>B i>j a>b 



2^ qaqBctj ■ Qi In 

a,B 



(4.15) 



We have slightly changed notation, compared to (4.14), notably the subscript m refers now 
to bions and e — to W^s. The interaction strengths are denoted by Kg for the electric charges' 
Coulomb interaction and ^ for the magnetic charges' interaction. The map between the 
UV-cutoff (a) values of Kg, ^m, Ve, Vm, and g, ^bion, iw, L, and T, is: 



Ve ^ (.WO^ , Vrn ^ CbionO^ , l^e{a) = tim{a) O 



27rLT 



(4.16) 



The electric charges (VF's) of type i (of charges qAoi^, qA = il) are located at R\ and the 
magnetic charges (bions) of type i (of charges qaQ^, qa = =tl) are located at i?^. Note that 
only e- and m- charge-neutral, with respect to [/(l)^"""^, configurations will give a finite 
contribution to the partition function. There are arbitrary numbers of electric (NJ,^) and 



-40- 



magnetic (NJ^^) charges of each type in the plasma, which are summed over in the grand 
canonical partition function. Since we are working at a center-symmetric point along the 
spatial circle, the fugacities of all Nc bions are equal (as well as those of the Nc W^-bosons). 
The dimensionless fugacities are defined as in (3.14) in terms of the bion and W fugacities. 
Similar to the fugacities of the different bions, the gauge couplings of all the unbroken C/(l)'s 
are equal. The "clock" symmetry of the extended Dynkin diagram responsible for the equal 
fugacities is preserved by renormalization and greatly simplifies the resulting RGEs. 

An important observation distinguishing Nc = 3 from the Nc > 3 theories is now due. 
From eqn. (4.8), it follows that for Nc = 3 the partition function (4.15) is invariant under the 
duality transformation: 



exchanging the fugacities of bions and VF's and appropriately inverting the coupling (note 
that for SU{3) Kg = t^m at any scale, see the next Section 4.2). As was already mentioned in 
the Introduction, this is because the interpretation of (4.8) is that different kinds of bions and 
W^s have identical interactions (which are only nearest-neighbor in the Dynkin diagram). We 
expect that the duality relation (4.17) is preserved by the RGE flow (and show that indeed, 
it is) and that its use should greatly simplify the study of the critical theory. 

On the other hand, since for Nc> 3 there is no relation between at ■ aj and Qi ■ Qj (see 
(4.2), (4.7)), the partition function is not invariant under electric-magnetic duality. Thus, we 
expect that the UV cutoff relation Kc = Hm is going to be violated by the renormalization 
flow. Thus, in addition to the RGEs for the two fugacities ye and j/m, we will have separate 
RGEs for Kc and Km- This expectation is, indeed, borne out by their explicit form given in 
the next section. 



4.2 The RGEs for the SU(Nc) magnetic bion/W-boson plasma 

Although there are different methods to derive the RGEs, we choose to follow the more 
intuitive approach introduced for the SU{2) case in Appendix A, and worked out for the 
general SU{Nc) model in Appendix C. 

Unlike the SU{2) self-dual model, the self-duality is lost in the SU (Nc) group with Nc > 3. 
Therefore, one does not expect the magnetic and electric coupling strengths to scale the same 
under the RG flow. This, in turn, forces us to define different strengths Ke and — (7^ — ) for 
the Coulomb interaction of electric and magnetic charges. The details of the procedure are 




12 



(4.17) 
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given in Appendix C. The final equations (C.6), (C.7), and (C.l), respectively, read: 



~db 



i=0,±l i=-l,0,l,2 



2 

p+i) 



j=0,±l,±2 j=-2 -1,0,1 

dye ( \ 

Notice the different number of terms in the sums over bion-pair and VF-boson pair contribu- 
tions to the (anti-) screening. This is because a given bion can interact with its next-to-nearest 
neighbors on the Dynkin diagram (i.e., it interacts with similar kinds bions as well as with 
four other kinds), while a given 1^-boson can only interact with the same kind of VF's and 
with its two nearest-neighbors. 

As already mentioned, A'c = 3 is an exceptional case. There, the sums over i in all terms 
in (4.18) are over 0, ±1 only. We can evaluate the r.h.s., using (4.8), (4.2), and (4.5), and 
cast the SU{3) RGEs in the form: 



dKe 

~db 



db \ Km. J \ k'L~^ 













2 J 




6 











ye, (4.19) 

ym- 



dye 

db 

dym 
db 

The UV relation Kg = is preserved by the RGEs (one can explicitly check that the second 
equation follows from the first, third and fourth). The RGEs are invariant under the e-m 
duality (4.17) {k^ — )• 12/Ke, ye -H- ym)- Further, we also observe that there is a fixed point 
for the Ke (or l/^e) equation at the e-m self-dual point Ke = = 2^/3, ye = ym- However, 
this is not a fixed point of the ye,m RGEs. Instead, both fugacities are relevant at the 
i^e,Tn = 2\/3 fixed point, at least according to the leading-order RGEs. While the existence of 
a self-dual point implies, as discussed in the Introduction, that Tc is the one corresponding 
to Ke,m — 2\/3, further study of the approach to criticality and the critical theory clearly 
requires extensions of the methods used here. 

We shall make very few comments on the SU (Nc > 3) case. The deconfinement transition 

2 

there also occurs in the region where both fugacities are relevant, in the window < Tc < 

2 

(equivalently 3 < Kg = K-m < 4 ) . It is also clear that studying the nature of the transition 
is more challenging than in the Ne = 2 case. This is because in the Nc > 2 case the RGEs 
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do not display a fixed point, at least to lowest order in fugacities, and because e-m duality, 
which is helpful to find the value of Tc even in models where small-fugacity RGEs break down 
(and even study the theory at Tc) is absent in QCD(adj). 

In the next section, we describe the construction of a lattice model dual to the Nc = 3 
Coulomb gas. The model is related to vectorial Coulomb gas models used to study the 
melting of 2d crystals with a triangular lattice [5] (not surprisingly, in view of the root lattice 
of SU{3)), by the addition of appropriate symmetry breaking fields. In this paper, we will 
only present the lattice model and will defer the study of its phase diagram for future work. 

4.3 The SU(3) Coulomb gas as a system of coupled lattice XY spins 

In this section, we map the Coulomb gas of bions and PV^-bosons for the SU{3) theory to 
the Z3 X Za-preserving "vector" spin model discussed in the Introduction, see eqn. (1.6). We 
provide two alternative, but equivalent, descriptions of the spin-system. 

4.3.1 First description 

The first description uses a two-component vector Ox = (0^, 6^), which is a periodic variable: 

4 = 4 + 27rai = 4 + 2^52 or 4 = 4 + 27rr^ , (4.20) 

whose periodicity is determined by (27r times) the root lattice Si G Tr of SU{3). The vector 0^ 
is associated with a site rr on a two-dimensional lattice. We shall take the partition function 
of the model to be defined as a path integral over the compact 9, Z = f T>6e~^^ (in the 
path-integration, 6 needs to be integrated over the unit cell of the root lattice T,.), where: 

3 3 
-PH= 5^^cos2z7i-(4+A-4) + X]Z]^^°^2(«^^-«^-i)-^^- (^-2^) 

x;fi=l,2 i=l X i=l 

Here, z7j are the three two-component weights of the defining representation of SU{3), Si — 
Si-i = Qi is the magnetic bion charge defined in (4.5). The factor of two in both terms 
ensures the periodicity (4.20) of 6?^ 

The continuum description of the lattice model (4.21) can be found by taking the naive 
continuum limit and keeping in mind that the two scalars are periodic according to (4.20). 

■^^In this footnote, we give an explicit basis for calculations and various properties of weights and roots 
of SU{3). As per our trace normalization, Tr{tatb) = |(5a6, the weights of the defining representation are 
i7i = (|, j:^), 1^2 = 2v^)' ~ roots are differences of weights and are given by 

ai = ui ~ U2 = (1, 0), a2 — ^2 — Us — (— |, and affine root as — us — ui — (— |, — Their length is, 
following normalization of the trace, normalized to one. The following relations are valid for general SU{N) 
and are useful in practical calculations. The weights, now represented by A'^ (A— l)-dimensional vectors Ui, 
obey: 

JV 

JV 
1=1 
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ill _ 

2 2iV' 



(4.22) 



Expanding the first term in (4.21) to quadratic order around the origin, '^'^ we obtain a 
canonical (up to an overall prefactor) kinetic term for the 9{x) field. Explicitly, denoting 
V^^x- = Ox+fi - Ox, we have: 



3 2 3 



x,ij, i=l x,^i a, 6=1 i=l 

^ Y.^y,Ox) ■ {V,Ox) -^-^jcPx {d,9f . (4.24) 



47r 

We shall take the partition function of the continuum version of the model to be defined as 
a path integral over the compact 0, Z = f T)0e~^^: 

-f3H = -^ [ {d^,ef + 2CY f cos2Qi-0 , = e + 2TrTr, (4.25) 

47r J^2 ^-^182 

where Qi = Si — ai-i is the bion charge. 

Expanding the interaction potential, exactly as in (3.22), we observe that the partition 
function of (4.25) precisely reproduces the magnetic bion interactions and their contribution 
to the partition function (4.15) with Q ~ ibion ~ Vm- Eqn. (4.25) also incorporates the W- 
boson interactions as well as the Aharonov-Bohm type e-m interactions, as in the discussion 
of the SU{2) case. The VF-bosons, once again, arise as vortices of the spin model. To see 
this, recall that in the continuum version is periodic by elements of the root lattice (4.20). 
Thus, a vortex at the origin can be written as: 

0{ip) = [n^ai + n^aa) ^ , (4.26) 

where n^jU? G Z and ip is the polar angle. Following Nelson [5], we take the fugacity of a 
vortex labeled by (n^,n^) at x to be given by:^^ 

y(n'{x)di)^ _ y{n\+n\-n^n2) _ (4.27) 

Thus, for small y, the leading-order contribution to the vortex partition function is due to 
vortices with winding number vectors denoted by ua, A = 1, 2, 3: 

ni = (1, 0), n2 = (0, 1), and = (-1, -1), (4.28) 



and they generate the weight lattice r„ of SU{N). The roots obey: 

di = Vi — Vi+i,Vi ■ oij =Vi- (vj — Vj+i) = \(5i 
di ■ dj = {vi — Vi+i) ■ {vj — i'j+i) = Sij — — |<5i,j-i , (4.23) 

and they generate the root lattice Vr- The root lattice is a sublattice of the weight lattice, and it is coarser. 
The quotient is T^/Vr = Zat. 

•^"Note that this is the unique minimum of the "kinetic" term in the fundamental domain of (4.20). 

''^We note that the lattice model produces a fugacity exactly of the form (4.27), with y determined by 
the details of the lattice cutoff; that this is so can be seen also from the continuum description, where the 
self-energy of a vortex is proportional to {n^dif' . In the Villain representation of the spin model, a fugacity 
term (4.27) can be added with arbitrary y. 
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and their corresponding anti- vortices. Note that the vortex with (ni,n2) = (—1,-1) corre- 
sponds to the affine root as = —di — 0.2 of the 5?7(3) algebra. It is pleasing to see that a 
T^-boson associated with the affine root — which is only present in the compactified gauge the- 
ory, and can be referred to as a Kaluza-Klein 14^-boson — arises naturally, with fugacity equal 
to that of the other two VK-bosons, from the spin-system as the vortex with charge assignment 
(— 1, —1). Each of the three vortices (4.28) gives the same contribution, proportional to y, to 
the vortex partition function. Any vortex with winding numbers different from (4.28) (and 
the corresponding antivortex) gives a contribution to the partition function proportional to 
a higher power of y and its effect is thus suppressed with respect to that of (4.28). 

Next, we can also easily see that the interactions of two vortices (4.26), separated by a 
distance r, with vorticities by n(0)*aj and n{ryaj, is given by: 

-K n{Oydi ■ ajn{ry ln{rT) , (4.29) 

which, when restricted to the leading vortices (4.28), reproduces the interactions between 
the three types of VK-bosons in (4.15). After also taking into account the e-m interactions, 
exactly as in the SU{2) discussion, we find that the affine spin-system (4.21) is equivalent to 
the gauge theory with partition function (4.15). 

The SU{3) /Tjs theory: The description (4.21) and its continuum form (4.25) are appropri- 
ate to the SU{3)/Z3 theory. They allow the introduction of dynamical monopoles of smaller 
charge than the bions, but do not allow the inclusion of particles of charge smaller than that of 
W^-bosons. This means that the theory does not possess a center symmetry, but a topological 
symmetry. The symmetry is Zg'^ x Z3 where the first factor is identified with Z3 discrete chiral 
symmetry and the latter is the Z3 topological symmetry associated with 7ri(S'J7(3)/Z3). As 
(4.21) shows, both Z3 symmetries are broken at low T, where the spin model is semi-classical. 
In the high-temperature phase, following direct analogy with the SU{2) discussion, none of 
the discrete symmetries is spontaneously broken. 

The SU{3) theory: In analogy with the SU{2) theory, we may construct a sigma model 
Lagrangian to describe the SU{'i) theory as well. In that case, the SU{2>) theory allows the 
introduction of dynamical matter fields whose charge is smaller than the VF-bosons', but does 
not allow the inclusion of monopoles whose charge is smaller than that of 't Hooft-Polyakov 
monopole. Consequently, the theory does not possess a topological symmetry, but just the Z3 
center symmetry. The discrete symmetry of this construction is "L^ x Z|. At low temperature, 
Zg'^ chiral symmetry is broken and Z| is unbroken. At high temperature, this is reversed: Zg 
is broken, and Zg"*^ chiral symmetry is unbroken. 

Electric-magnetic duals of (4.25,4.21) as well as a more detailed study of the phase dia- 
gram of the SU{2>) QCD(adj) theory will be considered elsewhere. 

4.3.2 Second description 

Here, we mention briefly a second lattice action which may be easier to simulate and which 
is equivalent to (4.21). Recall that the fields in (4.21) have a non-cartesian periodicity. 
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determined by 27TTr, the root lattice, and a canonical diagonal kinetic term for the associated 
continuum sigma model. It is also useful to provide a description in terms of two compact 
scalar fields on Mp that we denote = </>'(x) + 2tt {i = 1,2), which are both periodic by 

2tt. However, now the metric of the sigma model is non-diagonal. The lattice model is: 



/3H = — ^ cos(0i.+^ cos(0^.+^ -(1)1) + I sin((/.^+^ - (/)^.) sin((/)2_^^ _ (/.^ 

+ Y y (cos Scpl + COS 3(1)1 + cos 3(0;^ - c/)^)) , (4.30) 



where^^ y ~ Cbion ~ Um- The continuum version of (4.30) does not have a diagonal kinetic 
term. Instead, there is a metric gij = 1, 2) on the target space, encoding the interactions 
between vortices (i.e., the P^-bosons), gij = ai ■ aj. The continuum limit of (4.30) is: 

-fiH = -^ [ gijd^<P'd''^ + 2^bion I (cos 30i + cos 34!^ + cos 3(</.i - 0^)) . (4.3I) 

The Z3 X Z3 symmetry-breaking nature of the perturbation is also more manifest in this 
formulation. The charges (3,0), (0,3), and (3,-3) correspond precisely to magnetic bion 
charges in the frame. 

The relation of (4.30) to (4.21) (and of (4.25) to (4.31)) can be most easily elucidated by 
the change of variables: 

0i = + , 02 = -^1-^2/73. (4.32) 

The transformation (4.32) rectifies the fundamental domain (4.20) of 6 to that of (j)^''^ (0* = 
(j)'^ + 27r) and can be easily seen to map (4.31) to (4.25). The kinetic term of the lattice model 
(4.21) when expressed in terms of (p^''^, reads: 

_^^fc.n^ ^ £^cos(0i+^-</.i) + cos((/)2+^-</.2) + cos(0i+^-</.i-<Al+^ + (^2)J ^ (4_33) 

and differs from the kinetic term in (4.30) by an overall factor of 1/2 and the third term. 
However, we believe this difference to be inessential. First, it is easy to see that both (4.33) 
and the kinetic term in (4.30) have a unique minimum in the fundamental domain at V^cp^ = 
V^02 = and, hence, there are no "doublers" in either description. Second, the expansions 
of both (4.33) and the kinetic term in (4.30) around the minimum coincide to quadratic order 
in V ^(j)^'"^ , so that the models are identical in the Villain approximation; we expect that the 
difference (a four-derivative term) will not be essential in the continuum limit. 



5. Possible future studies 



Some possible directions related to the affine XY-spin model/gauge theory equivalence that 
are worthy of pursuit and would directly extend the results of this paper are listed below: 

^■^As promised in the Introduction, we note that the above action, with y = 0, is the one relevant in the 
study of melting of the two-dimensional triangular-lattice crystals [5] . 
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1. One can study the renormalization group equations for higher-rank QCD(adj) to next 
order in fugacities and search for fixed points, similar to studies of extended dual sine- 
Gordon theories [48-50]. 

2. In the SU{3) case, using 2d CFT methods, as in [4], one can focus on the self-dual point 
and attempt to infer the existence and nature of the fixed-point (critical) theory. 

3. A mean- field analysis could be useful to gain some idea of the behavior of the system 
in the region of large fugacities (to which the leading order RGEs appear to drive 
the theory). Recall that in the BKT transition, there is a tricritical point separating 
the small-fugacity continuous transition from a line of first-order transitions at large 
fugacity, revealed by a mean-field analysis, and observed in experiments and simulations, 
see, e.g., [53-56]. 

4. Numerical studies of the phase transition in the 2d lattice models with short-range 
interactions, dual to the long-range interaction Coulomb gas, should be possible. These 
would be useful to determine the order of the deconfinement transition. 

5. The study of non-equilibrium properties of thermal gauge theories is a subject of current 
interest, e.g., for RHIC physics. It would be of interest to extend the x §^ x 8\ 
framework of this paper to allow the study of the near-critical behavior of quantities 
such as the conductivity and viscosity. 

There are also other interesting directions, broadly related to the technology proposed here, 
which we intend to study more systematically: 

1. Orbifold/orientifold equivalences and spin systems for general QCD-like theories: In 
the large- limit of gauge theories, there arc exact orbifold and orientifold equivalences 
relating gauge theories with different microscopic matter content [58,59]. These equiv- 
alences have remnants in the small-iV semi-classical regimes, as emphasized in [25,60]. 
Therefore, by using double-trace deformations, our study should generalize in inter- 
esting ways to QCD with fundamental, anti-symmetric, and symmetric representation 
fermions. In particular, the semi-classical domains are also suited to check the univer- 
sality of the phase diagrams of QCD-like theories related through valid equivalences [61] 
satisfying the necessary and sufficient conditions [59] . 

2. Chiral gauge theory /spin system mapping: We certainly hope that the numerical cross- 
checks of the spin system/gauge theory equivalence will be successful and useful in 
vector-like gauge theories. If so, then it is worth generalizing the equivalence to chiral 
gauge theories, by using the techniques of [62,63]. 

3. Zero-mode Hamiltonian with fermions: In order to understand the meaning of the 

topological symmetry better, we have constructed the zero-mode Hamiltonian for the 
theory on x where S\ is parametrically smaller than T^. We aim to study the 
bosonic and fermionic zero-mode quantum mechanics for the quantum field theories 
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systematically. It seems to us that this framework may provide opportunities to relate 
the instantons with fractional topological charge proposed by 't Hooft [38] on toroidal 
compactification and the monopole-instantons or bions which appear semi-classically on 
xS\. The magnetic bions and monopoles, depending on the theory, provide reliable 
(i.e., semi-classically calculable) confinement mechanisms on x S^. On the other 
hand, ref. [64] provided numerical evidence for the relevance of objects with fractional 
topological charge on the torus, for which analytic solutions are not known. We aim to 
gain a better understanding of their interconnection. This project is ongoing. 

4. Zero-mode Hamiltonian and 9-angle dependence: We expect that our framework on M x 
X 'B\ will provide new insights into the topological 0-angle dependence of observables 
in gauge theories. This project is ongoing. 
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A. RGEs for the SU(2)(adj) magnetic-bion/W-boson Coulomb plasma 

Consider the partition function of a form more general than (3.12), namely that of a neutral 
gas of unit "electric" and charge-p "magnetic" charges, where the interaction strengths are 
given by k for the electric charges Coulomb interaction and ^ for the magnetic charges; in 
addition, there is the appropriate Aharonov-Bohm phase (3.11). The dimensionless fugacities 
are ye and ym, respectively, and a is the "lattice spacing". The partition function reads: 



Z EM,p 

(y^^2N^ (y,\2Ne 2N^ 2Ne 



(A.l) 



X exp 



- QaQb In — — + K qAQB In ^— — ^ + qBqa&(.RB - Ra) 

a>b A>B a,B 



The map of (A.l) to the bion/W-gas partition function (3.12) is straightforward, using the 
change of variables (3.14). Our goal is to derive the RGEs (3.16) to leading order in fugacities 
ye, Vm, describing the behavior of the three parameters k, ye, ym under course graining of the 
partition function. These equations are well-known and we only present our derivation to 
make our exposition self-contained (we believe ours is the simplest derivation to leading order 
in fugacities; for higher-order calculations, see [48]). 

Consider first the renormalization of the fugacities. They can be most easily obtained 
by considering the lowest-order contributions to Zem,p- For example, the {Nm,Ne) = (1,0) 
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term reads: 



ym{a) 





ri - r2 




a 



(A.2) 



where we have indicated that the fugacity on the l.h.s. is the one appropriate for lattice 
spacing a, the ri^2 integrals were performed from a to the size of the system L. Demanding 
RG invariance of the contribution (A.2) to the partition function: 



L 



(A.3) 



and differentiating (A.3) w.r.t. 6 at 6 = 0, we obtain: 

-,2 



dym 
db 



(A.4) 



the third of eqns. (3.16) given in the text (with ym = y; note also that to leading order in the 
fugacity we can neglect the change of k in (A.3)). The RGE for Ue, the second of eqns. (3.16), 
follows by duality or from similarly considering the {Nm,Ne) = (0,1) contribution to the 
partition function: 

^ = (2-f).„. (A.5) 

It is easy to check that the scaling of the fugacities (A.4) and (A.5) also ensures that the 
entire partition function is invariant under course-graining. 

To find the effect of screening and anti-screening on the interaction between magnetic 
or electric charges, consider for definiteness the interaction between two external magnetic 
charges with q = +1 at x and q = —1 at y. To leading order, ignoring any fluctuation 
contribution, the interaction between external magnetic charges is given by: 



-7-T- In \x—y\ 



(A.6) 



where we have indicated that the strength k is the tree-level value appropriate for a cutoff 
a. Our goal now is to compute the effective strength of the interaction appropriate to a 
course grained system (i.e., find K{e''a)). In other words, we want to find the change of the 
interaction (A.6) due to fluctuating electric and magnetic dipoles with sizes between a and 
e^a. The expression for their effect takes the following form (where we use the shorthand 
notation riy = |ri — y\ and @x2 = Q{x — 7^2), etc.): 



^Vm'^'^ I ^2^^^2^^g-7^(lnr:,a+lnri2+lnr:,i-lnr^2-lnriy+lnr-2a) 



+ 



1 + 



Vm 



d rid r2e «(°) 



+ 
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The meaning of the various terms above is as follows: the term ~ y,^ represents the contribu- 
tions of fluctuating magnetic dipoles (a positive charge at r2 and a negative charge at ri) to 
the interaction between the external magnetic charges and the term ~ yl similarly represents 
the contribution of an electric dipole (a positive charge at r2 and a negative charge at ri). 
In each case the energy is equal to the sum of the energies of interaction of all charges in the 
configuration, and to the order of the calculation involves six pairwise interactions in each 
case, as shown explicitly in (A. 7). The last term in (A. 7) is the \/Z normalization factor of 
the expectation value of the potential between the external charges, also written to leading 
order in fugacities. When this term is expanded to order ylyy^, the expression for the effective 
interaction takes the simpler form: 



•s.(e'>a) 



- ^ In ra: 



1 + 



Um 



' ^ (In r^,! -In ra;2 -In riy +ln r^y) 



K In ri2 



-ip(Qxl-&x2-Qly+Q2y) 



1 



(A.8) 



which, as appropriate, only takes into account the interactions between the dipoles and the 
external charges (the coupling k on the r.h.s. is the one appropriate at the scale a, but 
beginning with (A.8), we do not explicitly indicate this for brevity; also for the time being 
we do not explicitly show the limits on the integration over the positions of the fluctuating 
dipoles ri^2)- 

Recall now that we are interested in computing the effect of magnetic and electric fluctu- 
ating dipoles of small sizes on the effective interaction between the external charges. To this 
end, introduce centre of mass (R) and relative (r) coordinates: 



\nrxy 



+ 
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[1 + 



d^re 
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-2^1nr 



(A.9) 
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;+lnr 



R-y, 



d^R 



-ip{&R- 



-e 



R-y 



r+e 



R-y,- 



and note that for small r, f{rR_^^r_) - /(r/j„a, _ | ) - /{rji^y^r) + f{rR_y^^^_ ) ~ f • VnifirR^^) - 
f{rji^y)) (we have in mind that the function /{vr^x) stands for either ln{rji^x) or Qr,x)- Next, 
we note that expanding the exponents and averaging over the direction of r, the contribution 
linear in r vanishes, while the second order contribution becomes: 



-^^Inr^ 
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Above, we explicitly indicated that we are only taking into account the (anti-)screening 
contribution of electric and magnetic dipoles with sizes between the old and new cutoffs a 
and e^'a. Now we use the identity {V r In rji^x)"^ = (VrQr^x)'^ — an expression of the "electric- 
magnetic" vortex-charge duality in 2d — to recast (A. 10) as: 



rjTT ^^^xy 



P'^ 1 



(A.ll) 



+ 



2k? 



2/e \ 2 vrp^ 



drr e 



3 K In r 



Finally, we integrate I{x,y) = J cPR R{\nrRx — ^'^i'R,y)j by parts, using the neutrality 
of the system to discard the boundary term, to obtain I{x,y) = Airlnrx^y — 47rlnro.o, and 
ignoring the constant term, which does not contribute to the potential between the external 
charges, substitute back into (A.ll) to obtain: 



e " 



" [l-p"^ In 



(A.12) 



V 



' drr'^e - 



Re-exponentiating, we find that the effect of electric and magnetic fluctuating dipoles of sizes 
between a and e^a renormalize the strength of the interaction between external magnetic 
charges as follows: 
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Ve 



27r^ / drr^e 



3„— Klnr 



y^\ 2 27rV 



/ ^^(a) 



. (A.13) 



Finally, we differentiate w.r.t. b and take 6 = 1 and, remembering that the interaction is 
really In - , obtain: 

(A.14) 



d f^] o 2 / 2 2 P 



or equivalently. 



o 2 /' 2 2 2 2\ 

K = 27r [p y^- K 
For p = 4, this is exactly the first of eqns. (3.16). 



(A.15) 



B. The critical point of SU(2)(adj) and the approach to Tc 

We begin by rewriting the RGEs (3.16) of the Z4 model in terms of the variables: 

X = Y = A-Ky, Z = Any , 



(B.l) 
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keeping up to quadratic terms in X, Y, Z, which is sufficient for studying the approach to the 
fixed point: 



Y = -XY , (B.2) 
Z = XZ . 



These equations have two integrals of motion: 

YZ = c , (B.3) 

y2 + - = c , 

hence the RGEs can be solved in quadratures. 

See Figure 5 (Section 3.2) for an illustration of the RG flow (B.2), which we will now 
analyze in some detail. 

B.l The critical RG trajectory and determination of Tc 

We begin by considering trajectories that end on the critical line. To this end note that c' is 
determined by the initial (UV cutoff scale) values of the fugacities, c' = Y{0)Z{0) = IGir'^yoyo, 
with yo,yo given in (3.14). The fixed point of (B.2), whose location on the fixed line n = 
4,y = y is determined by the initial conditions, is at X = 0, y = = -v/c?- The last equation 
in (B.3) then implies that the fixed point trajectory is the one for which c = 2c'. We can now 
solve (B.3) for X, using this value for c: 



X ^ , ,B,4, 

We can now use this to integrate the equation for Y{t) from (B.2), taking (without loss of 
generality) 1^(0) > ^/d and a positive sign for X: 

Y{t) = Vdianh (^Vc't + tanh"^ ' ^^'^^ 

showing clearly that this trajectory approaches Y = \fd as t — t- oo. The corresponding 
solution for Z is: 



tanhf\/?t + tanh^i 

The solution for X{t) is given by (B.5) substituted in (B.4). As i — ?• oo this clearly approaches 
the X = fixed point, while as t — )• 0, we have: 



Note that l^(i) and Z{i) monotonically change from 1^(0), Z{^\ respectively, to y d = Y = Z 
as t changes from to oo. If the initial values Z(0) and ^(0) are small enough, and close 
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to the fixed-point value \fd ^ ^(i) is also small for all t G (0,oo), so that the approximations 
made in arriving at (B.2) hold. Thus, the entire RG trajectory to the critical point lies within 
the validity of the perturbative RGEs (B.2). 

We conclude that the above analysis of the approach to the critical point is reliable, 
provided the initial values of the fugacities and «; — 4 are small enough — as they are for our 
small-L S'C/(2)(adj) theory. The solutions to the RGEs (B.5, B.7) ending at the critical point 
are then reliable all along the flow — in contrast with the flow in the Z2 model (and more 
generally, in the Zp<4 models) relevant, e.g., in the deformed- Yang-Mills study of [12]. 

Equation (B.7) is a constraint on the initial conditions for a trajectory to end at the 
critical point. Thus, it can be used to determine Tc- Recalling the definitions (B.l), (3.14), 
as well as (3.8) and (3.6) for the VF-boson and bion fugacities, using miy = ^, and taking 
the lattice spacing to be the bion size (3.2), we have: 

^ - 4 ^ ..rl ((2„, . 4 - .-¥^<'-') . (B,8, 



where g = g{L) is the four-dimensional running coupling at the scale L, c = y -^3 — , and A is 
the unknown numerical coefficient of the bion fugacity (only the dependence on the coupling 
is known, see [35] for details). 

The point of equation (B.8) is that the critical temperature can in principle be precisely 

calculated by careful matching of the effective 2d theory to the 4d UV theory in the small-L 

2 

regime. Clearly, one can iteratively solve (B.8) to find the deviation of Tc from (equiva- 

lently, the deviation of kq from 4). However, since the fugacities are exponentially small, it is 

2 

clear that Tc = £^ is an exponentially accurate estimate — the shift of Tc is easily seen to be 

f(\ o7y 

of order ^-j^e where f{g) has a known leading-order dependence on g but is itself only 
known up to a numerical coefficient. 



B.2 The correlation length as T — > Tc 

Let us now analyze the behavior of almost critical trajectories. For the critical trajectory, 
we have from (B.7) that X(0)c = ^(0)c — ■Z^(0)c. Now, consider a trajectory with X(0) = 
X{Q)c — A and use the definitions X = k = ^^^^ , and Tc ~ to find: 

T — T 

A = -4—^. (B.9) 

c 

A trajectory with the chosen value of T has to also obey (B.3), in particular, at t = 0: 

c = Y{Qf + Z(0)2 - X(0)2 ~ Y{Q)l + Z{<d)l - ^(0)c + GA = 2c + GA, (B.IO) 

where we used c = 2c' for the critical trajectory, expanded for T near Tc (small A), and 
denoted by G the coefficient of the first term in the expansion of Y"^ + — X"^ near Tc- The 
value of G will not play a role in determining the critical properties (it can be calculated using 
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(B.8)). Let us now study the non-critical RG trajectory of (B.2), for which Y{t) is given by: 



Y{t) 

dY 

-t. B.ll 



y(o) 

For c given by (B.IO) we rewrite (B.ll) as: 

Y{t) 



^y4 _ c y2 + 



/^2 



'^^ -t . (B.12) 



y(o) 



^(y2_ c')2-GAy2 



It is possible to integrate (B.12) explicitly, but for the purposes of studying the divergence 
of the correlation length as A — )• 0, the argument below is sufficient. Clearly, for A = 0, the 
trajectory is given by the critical solution (B.5). It is clear that for small A, the trajectory will 
be essentially the same as the critical one up to values of RG "time" when \Y{t^)'^ — c'p ~ 
|GA|y(i*)2. Thus, we can use the critical solution up to times of order t*, which we can 
estimate by requiring that the two quantities below are of the same order of magnitude: 

\GA\Y{uf ~ |A|tanh2(^/?t,)~|A|(l-2e-2^**) , 
\Y{uf - cf ~ (1 - tanh2(V?t,))2 ^ g-^v^** , (B.13) 

where we assumed large t* and neglected numerical coefficients. Equating the two lines in 

1 

(B.13), we have ~ — ln(| A| ), hence the correlation length C, ~ e** scales as A ^ as: 

^ ~ \A\~T77 , (B.14) 

i.e., with a non-universal exponent, whose value depends on the fugacities of the bions and 
W^-bosons, y and y. The approach to Tc is governed by their values at the fixed line (equal to 
\/~d), for which, according to (B.3) we have c' = IGvr^T/oyo (where the subscript denotes the 
UV values of the fugacities). Thus, we can write, using the map (B.l): 



C~ \T-Tc\ i«-v^ , (B.15) 

thus we have v = ——p. — \ — ^ , in the usual convention C ^ \T — Tc\~^ . Recall that the theory 
is in a gapped phase on both sides of Tc (related by the electric-magnetic duality) and hence 
the scaling holds upon approaching Tc from both sides. The singular part of the free energy 
is: 

1 

F ~ ~ |r - Tcl^'v^^ . (B.16) 

Thus, for small fugacities {yo ~ e in the small-L S'[/(2)(adj) case), the transition in the 
Z4 gas is of quite high (but finite, unlike the Kosterlitz-Thouless transition) order determined 
by the integer part of 2v. 
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C. RGEs for the SU(Ne) case 



In this section we consider the RGEs in the SU{Nc) case. At temperatures much smaller than 
M\Y and larger than the photon mass, the system can be describes as a 2D neutral gas of Nc 
distinguishable copies of "electric" and "magnetic" charges. To derive the RGEs to leading 
order in fugacities ye and ym, we now closely follow Appendix A. The renormalization of 
the fugacities can be obtained from the lowest-order contributions of the partition function. 
Taking the contribution to (4.15) with N^_^ = = 1 (for some fixed p), we obtain: 

dye _ ( ay -dp \ _( 

while taking Nl^_^ = = 1, we obtain: 

dym _(^_2Q^\^^/_^ forAr,>3, (C.l) 



db \ 

where we indicated that the last equality in the second equation holds for iVc > 3 only. 

To derive the RGE for Kg, we consider two external electric charges of the type p, where 
1 < p < with charges q = +1 at x and q = —1 aX y. Ignoring the fluctuating contribution, 
we find that the interaction between these two charges is given by: 

^-ap-dpKe{a)\n\x-y\ ^ (C.2) 

Now, we take into account the contribution from the fluctuations of the electric and magnetic 
dipoles of the various types. The expression for the net effect takes the form: 

j (prid'^r2 ^Ke(a)dp-dp+i[lnrj;i-\nrx2-inryi+\nry^)-Ke(a)dp+i-dp+ilnri2 
^ i={0,±l} 

^ (^ymY f (frid'^r2 ^ ^'iiap-Qp+,{e^i-0^2-Qyl+&v2)-jr^Qp+i-Qp+ilnri2 
^ j={-l,0,l,2} 

where we used the fact that while every W boson interacts with itself and its two nearest- 
neighbors in the Dynkin diagram, a W^-boson interacts with four kinds of bions, as per 
(4.7) and (4.2). The partition function Z to leading order in fugacities reads (omitting the 
contributions of particles that do not interact with the external ones as they do not contribute 
to the (anti-)screening): 

" i={0,±l} 
" j={-l,0,l,2} 
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Expanding the Z ^ terms to order ^ we find: 



i={0,±l} 

+ (J^y J (fri(fr2 e^^^'^''+'''^^+'^'""'^ l^^2idp-Qp+,{e,i~ex2^eyi+ey2) _ i 



i={-l,0,l,2} 



Using the center of mass R and relative r coordinates and proceeding as we did in the RGEs 
of SU(2) we obtain: 



^-dp-dpKe{e''a)lnrxy _ ^-dp-dpKelnr^y l^]^ -)- 47]-lni 



xy 



TTK.l{a) /ye ^2 



2 \a 



_r, 1 1 J a 



^j,^3^-Ke{a)dp-dp Inr 



i=0,±l 



i=-l,0,l,2 



(C.4) 



Re-exponentiating we find: 



,2/ N fV' 



i=0,±l ■^'^ 



^j,^3^-Ke{a)dp-dp Inr 



i=-l,0,l,2 

Differentiating w.r.t. h we obtain: 



o-Qn Inr 



(C.5) 



db 



(C.6) 



j=0,±l 



i=-l ,0,1,2 



One can redo the same analysis for two external magnetic charges of strengths Qp and —Qp 
located at x and y to find: 



d ( 1 



4y2 



db\KmJ Qp-Qp 



j=0,±l,±2 



i=-2, -1,0,1 



Equations (C.l), (C.6), and (C.7) constitute the RGEs for the SU{Nc) case with iVc > 5 (the 
restriction comes from having to have five bions that interact with each other; for Nc = 3,4 
the equations will require a minor modification). 
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